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We present evidence that tunneling processes in near-integrable systems are enhanced due to 
the manifestation of nonlinear resonances and their respective island chains in phase space. A 
semiclassical description of this "resonance-assisted" mechanism is given, which is based on a local 
perturbative description of the dynamics in the vicinity of the resonances. As underlying picture, 
we obtain that the quantum state is coupled, via a succession of classically forbidden transitions 
across nonlinear resonances, to high excitations within the well, from where tunneling occurs with 
a rather large rate. The connection between this description and the complex classical structure of 
the underlying integrable dynamics is furthermore studied, giving ground to the general coherence 
of the description as well as guidelines for the identification of the dominant tunneling paths. The 
validity of this mechanism is demonstrated within the kicked Harper model, where good agreement 
between quantum and semiclassical (resonance-assisted) tunneling rates is found. 
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I. INTRODUCTION 



Since the early days of quantum mechanics, tunneling has been recognized as one of the hallmarks of the wave 
character of microscopic physics. The possibility of a quantum particle to penetrate an energetic barrier - or, more 
generally, a dynamical barrier in phase space, as termed by "dynamical tunneling" pj - represents certainly one of 
the most spectacular implications of quantum theory and has lead to various applications in atomic and molecular 
physics as well as in mesoscopic science (for instance in the context of spintronics devices (^] to mention just one recent 
example). Despite its genuinely quantal nature, however, tunneling is strongly influenced, if not entirely governed by 
the structure of the underlying classical phase space. Tuning the classical dynamics from an integrable to a chaotic 
one substantially modifies the tunnel coupling - not only if the chaos affects the classically allowed motion (a case 
that was studied in detail in (^]), but also if the destruction of invariant tori is restricted to phase space domains deep 
inside the classically forbidden regions. 

Tunneling processes take place in a variety of contexts. One may for instance be interested in the coupling of a 
metastable state to a continuum, as, e.g., in the case of the a- particle decay. The tunneling process manifests then 
in form of a finite life time of the state, or equivalently, in a finite width of the corresponding spectral resonance. 
Another situation concerns systems that possess a discrete symmetry P, in such a way that Rabi oscillations can take 
place between states that are localized in symmetric regions of the phase space. The effectiveness of tunneling is then 
related to the period of the Rabi oscillations, or equivalently to the splittings between the corresponding symmetric 
and antisymmetric eigenmodes. 

Here, for the sake of definitcness, we shall restrict ourself to this latter configuration, though our study could be 
applied with minor modification to other contexts such as the computation of resonance widths. To be more specific, 
we shall consider systems the classical dynamics of which exhibits, due to the discrete symmetry, two congruent but 
separate regions of regular, bound motion in phase space. Semiclassical Einstein-Brillouin-Keller (EBK) quantization 
then provides a prescription how to construct local eigenmodes on quantized tori within each of the regions. These 
"quasimodes" fulfill the Schrodinger equation, possibly to all powers in ft, but do not take into account the behavior of 
the wavefunction in the forbidden domain, and are therefore characterized by equal energies for any pair of symmetry 
related quantized tori. It is the nonclassical coupling between the quasimodes that lifts this degeneracy and selects 
the symmetric and the antisymmetric linear combination of these states as true eigenstates of the quantum system. 
Provided the classical dynamics is integrable, the tunnel coupling that produces this level splitting can be reproduced 
by the analytic continuation of the torus structure into complex domain. In case of a one-dimensional motion 
generated by the usual p 2 /2 + V(x) Hamiltonian, this amounts to calculating an imaginary action integral of the type 
/ y/2(V(x) — E)dx in between the two wells. 

Such a dynamical tunneling process also takes place in nonintegrable systems, but with substantially different 
effectiveness. By applying a nonintegrable perturbation on the dynamics in such a way that an appreciable chaotic 
layer is introduced in between the two regular regions, the tunnel coupling between the quasimodes is significantly 
enhanced with respect to the integrable system M~Q- Moreover, the tunneling rates do no longer exhibit a well-defined 
scaling with Planck's constant (which would be an exponential decrease with 1/ft in an integrable system), but may 
undergo huge quasierratic fluctuations for small variations of ft [Q-|6| . 

These phenomena are traced back to the specific role that chaotic states play in such systems fsj-fTj In contrast 
to integrable dynamics, the tunnel doublets of the localized quasimodes are in a mixed system no longer isolated in 
the spectrum, but resonantly interact with states that are associated with the chaotic part of phase space. Due to 
their delocalized nature, such chaotic states typically exhibit a considerable overlap with the boundary regions of 
both regular wells. They may therefore provide an efficient coupling mechanism between the quasimodes - which 
becomes particularly effective at values of ft (or any other external parameter) at which one of the chaotic level is 
shifted exactly on resonance with the tunnel doublet. This interpretation was basically confirmed by a statistical 
analysis of the fluctuations of the level splittings, which were shown to be in good agreement with the corresponding 
prediction that a random matrix model of the chaotic part of phase space would yield fj|J§. Evidence for chaos- 
assisted tunneling was found in theoretical calculations of the life times of nondispersive wave packets in microwave 
driven hydrogen || and, most recently, also in experiments on ultracold atoms that were stored within periodically 
modulated optical lattices ||io|-|l2||. Also for electromagnetic systems, in particular within microwave billiards {p| and 
optical microcavities [ fl4| , it was verified that (ray optical) chaos has a strong impact on wave tunneling. 

Apart from this qualitative insight and its quantitative statistical modeling, however, little is known about chaos- 
assisted tunneling from the semiclassical point of view. In particular, there exists no practicable semiclassical tool 
up to now for estimating the tunnel coupling from a torus within a regular island to the chaotic sea (which would 
provide valuable information about the magnitude around which the tunnel rates fluctuate) . This is essentially due to 
the fact that in nonintegrable systems invariant tori are generally characterized by a restricted domain of analyticity 
and can therefore not be continued far enough into the complex domain to reproduce the tunnel coupling |Iq,|l6|] . As 
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has been shown by Shudo and coworkers (T^|, it is in principle possible to semiclassically reconstruct the tunneling 
tail of the time evolution of some initial wavefunctions within a mixed system by summing up the contributions of 
all complex classical trajectories that fulfill appropriate boundary conditions. This requires the analysis of a highly 
complicated manifold of initial conditions in complex phase space, where the Stokes phenomenon plays a crucial role 
for the selection of contributing trajectories [ p~7[ . Though successful from the quantitative point of view, this method is 
limited to initial values problems, and it is not evident how much insight it may provide for an intuitive understanding 
of the properties of cigenfunctions. 

The present contribution can in some sense be seen as a complementary approach with respect to the work of 
Shudo et al. jl7|]. We shall, however, not tackle directly mixed regular-chaotic dynamics, but consider a conceptually 
simpler situation, namely the nearly integrable regime in which chaos is not yet fully developed. In this context, 
it has been pointed out by Bonci and co-workers p8| that classical resonances should play a crucial role in such 
systems. As nonlinear resonances between two classical eigenmodes arise whenever their frequencies are related by 
a rational multiple, they generally lead to quasi-degeneracies of the associated levels in the spectrum. This may 
therefore induce significant couplings to higher states within the regular region, from which tunneling occurs with a 
much faster rate than from the original state. Besides the aspect of near-degeneracies - which alone is not sufficient 
for a quantitative understanding of the role of resonances, as level crossings arise also in integrable systems without 
introducing any modification - it is also important to take into account the magnitude of the effective matrix element 
that couples the near-degenerate states. This quantity is directly related to the structure of the underlying classical 
phase space. Within the action-angle variable representation of the unperturbed integrable dynamics, a nonlinear 
resonance generally manifests in form of a pendulum-like perturbation, and the associated modification of the torus 
structure in its vicinity is then responsible for couplings across the resonance, the strength of the coupling being 
determined by the extension of the pendulum structure in phase space ]l9|-p2[ . 

Our previous publication p3| has shown that the combination of these two aspects, namely the near-degeneracy and 
the coupling induced by the resonance, form the basis of a mechanism that governs tunneling in the near-integrable 
limit. This study was performed within the "kicked Harper" model, a one dimensional time periodic system governed 
by the Hamiltonian 

oo 

H(p, q, t) = cosp + t 5(t — nr) cos q (1) 

n— — oo 

the dynamics of which is equivalent to the symplectic map 

p' = p + t sin q 

q' = q — t sinp' (2) 

which describes how the phase space variables (p, q) evolve from time t = nr (more precisely, from the time immediately 
before the kick) to time t = in + 1)t. This model has proven its usefulness in the context of many different aspects 
related to quantum chaos [p4|-p6| (including also dynamical tunneling p^|). Our study was restricted to a relatively 
small value r = 1 of the perturbation parameter, for which the classical dynamics in nearly integrable. The quantum 
tunneling rates that are obtained at this perturbation strength are shown in FigJ^. We see that, despite a seemingly 
"regular" phase space (shown in Fig. |J), they are nontrivial and exhibit similar features as in the case of a truly 
mixed regular-chaotic system: Even for rather small deviations from integrability, the tunneling rates may, in the 
semiclassical regime, become appreciably enhanced with respect to the integrable limit (by a factor that may reach up 
to ten orders of magnitude in the case that we have considered in |23|]) and do not follow a monotonous exponential 
scaling with 1/h. 

As key ingredient to understand such a behavior, we have introduced an integrable approximation of the kicked 
Harper map (in an analogous way as in p8||) which, even in the case of moderate perturbations, provides a reasonable 
description of the nearly integrable motion on the invariant tori. Expanding the kicked Harper eigenfunctions within 
the eigenbasis of this integrable approximation allowed to unambiguously identify resonances as the source of modi- 
fications in the tunneling tail of the eigenfunctions. A quantitative reproduction of the tunneling rates in the kicked 
Harper, the accuracy of which is visible on Fig. [l], was then achieved through a quantum perturbative treatment of a 
local effective Hamiltonian, which is formally derived via secular perturbation theory of the classical motion |29|| and 



was in practice obtained via the Fourier analysis of the separatrix structure associated with the resonance 23 1. 

The combination of these "tools" has evidently proven successful for the identification of the underlying mechanism 
as well as for a low-cost (with respect to computer memory) calculation of tunnel splittings that would otherwise 
be accessible only through a full quantum treatment of the problem. However, the justification of the resonance- 
assisted mechanism presented in j23fl was mainly based on the demonstration of its quantitative predictive power 
for reasonably small values of h within the kicked Harper system. In this paper, we would like to go further in the 
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understanding of the tunneling process in the nearly integrable regime. A central question that we shall address is to 
which extent the resonance-assisted tunneling mechanism we propose should in general be the dominant one, and what 
modifications are to be expected as the system is pushed deep in the semiclassical regime. A general issue that underlies 
these interrogations is the fact that the approach we propose is based on a combination of perturbative techniques 
(both classical and quantum) and semiclassical concepts, and therefore involves essentially two small parameters: 
the perturbation strength which enters in a purely algebraic way into the coupling terms, and the quantum coarse 
graining h on which these terms depend both algebraically and exponentially. Although we obviously do not intend 
to attain anything like mathematical rigor, our goal in this paper is to give evidence that the global picture that 
underlies resonance-assisted tunneling "makes sense" and, on a more practical tone, may lead to guiding rules for the 
identification of dominating terms in the tunneling mechanism. 

To reach this objective, we shall see that it is useful to provide a more geometric vision of resonance-assisted 
tunneling. This means on the one hand that we shall emphasize the connection between the coefficients that describe 
the strength of the coupling and the complex structure of the underlying integrable approximation. On the other hand, 
we shall see how the coupling via a nonlinear resonance can be considered as dynamical tunneling process, in very 
much the same spirit as the coupling between quasimodes on symmetry-related invariant tori. However, the effective 
topology of complex tori that the quantum system encounters in order to undergo the tunneling transition sensitively 
depends on the quantum coarse graining. For rather large h, a direct connection between the quasidegenerate tori of 
the two wells is "seen" by the quantum system. Deeper in the semiclassical regime, the tori rather appear as being 
connected, via one or several resonances, to higher excitations within the well, from where a transition across the 
separatrix is associated with a rather low imaginary action. 

Our study will be restricted to one-degree-of-freedom systems subject to a time-periodic perturbation with period 
r or frequency to — 2-k/t. We denote by H(p,q,t) = H(p,q,t + r) the quantum Hamiltonian and by H(p,q,t) its 
classical limit. The classical phase is most conveniently visualized by means of a Poincare surface of section in time 
domain - i.e., by the area-preserving map 

T ■ (Pu q t ) ' — ► (Pt+r, qt+ T ) (3) 

that describes the evolution of the phase space variables (p, q) from time t to time t + r. Quantum mechanically we 
shall, in analogy, consider the quantum propagator 

U = exp(j_^ H(t)d?j , (4) 

and study its eigenfunctions \ipk) and eigenphases 4>k, defined by 

Ufa) = e*>\il> k ) . (5) 

Whenever an illustrative example appears appropriate, we shall make use of the kicked Harper Hamiltonian |24J] Eq. (|l|) 
the Poincare map of which is given by Eq. (||). We shall, however, try to keep the discussion as general as possible 
in order to allow an application also to other time-periodic tunneling problems such as the driven double well p|]3"c|l 
or the effective Hamiltonian |l2| ] that was employed in the context of the recent dynamical tunneling experiments in 
cold atoms 0,0. 

To lay firm foundations, we begin in Sec. with a brief review of what we like to name "regular tunneling" - i.e., the 
attempt to semiclassically describe tunneling by a direct analytic continuation of the invariant tori into the complex 
phase space. We shall argue, however, that this concept is, strictly speaking, limited to exactly integrable systems 
and breaks down when a small nonintegrable perturbation is applied. This naturally leads to the qu estion of how 



nonlinear resonances influence tunneling, which we shall discuss in Sec. III. We shall begin, in Sec. HI A, with a formal 



description of the clas sical dynamics in the vicinity of a nonlinear resonance, based on secular perturbation theory, 



and use then, in Sec. [II B , quantum perturbation theory as well as semiclassical WKB theory to study transitions 
across the resonance. The practical calculation of the coupling co efficien ts that parametrize this description, and a 
discussion of the ge neral properties of their scaling, is given in Sec. [II Cj . Plugging these basic elements together, we 



then obtain, in Sec. [II D, a satisfactory semiclassical picture of how tunneling proceeds in presence of one or several 



resonances at given value of the quantum coarse graining. To demonstrate its feasibility as well as to verify basic 



assumptions that have been made in the course of its derivation, we finally return, in Sec. IV, to the particular case 
of the kicked Harper Hamiltonian, in a parameter regime where its classical dynamics is nearly integrable. 
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II. "REGULAR" TUNNELING 

A. Tunneling in integrable systems 

For one-dimensional time-periodic systems, integrability can be defined by the existence of a function J-(p, q) that 
is conserved by the Poincare map T describing the evolution of (p, q) from time t to time t + r. This can be shown 
to be equivalent to the existence of a r-periodic canonical transformation (p,q) i— ► (p(p,q,t),q(jp,q,t)) such that 
the Hamiltonian in the new coordinates is time independent f3lf| - in which case the conserved quantity is simply 
the energy. Without loss of generality, therefore, we discuss in this subsection the properties of time-independent 
Hamiltonians H = H(p, q). 

Integrability quite naturally yields a great number of simplifications. Due to the existence of a constant of motion, 
the iterates by the Poincare map of a given point in phase space lie on an invariant curve T (see, e.g, Fig. ||) which 
we call, in analogy to higher dimensional systems, a "torus" throughout this paper. It will be convenient to use the 
action-angle variables (/, 9) associated with H. For a given phase space point (p, q) on the invariant torus T, the 
action is defined by 

1 = f P d 1 ( 6 ) 



2?r Jr 

and corresponds, up to the factor 2ir, to the area that is enclosed by the torus in phase space. The angle 8 £ [0, 2ir[ 
represents the conjugate variable and corresponds to the propagation time that elapses from a given reference point 
on r up to the point (p, q) (normalized in such a way that 9 = 2n after one full round-trip). Expressed in these new 
variables, the Hamiltonian is, by construction, a function of the action only: 

H(I,9)=H(I). (7) 

(In order not to overload the notation, we shall use the same symbol H for the Hamiltonian in the original phase 
space variables (p, q) and in the action-angle variables (1,8)). 

Quantum mechanically, the time-invariance of the Hamiltonian H implies that the propagator U of the wavefunction 
from time t to time t + r (Eq. |]) is simply given by U — exp(— i_Hr/S). Its cigenfunctions tpk are then also the 
eigenf unctions of the Hamiltonian, and the associated eigenphases are related via (j>k = Ekr/h to the eigenenergies 
Ek of H. They can, moreover, be semiclassically constructed using standard EBK theory. More precisely, the 
semiclassical wavefunction that is associated with an invariant curve T is defined by 

* m w= _ e ' r Jyy^ + ,' r M) (8> 



\ 



I J'o 

r 



where v(xo,x a ) is the [algebraic] number of vertical tangents that are encountered by T between the phase space 
points xq and x a [ p2| . ip(l) can be properly defined (i.e. is mono- valued) if and only if the action enclosed by the 
curve r fulfills the quantization condition 

I = hj >pdq = h i K k+ ^) =Ik (9) 

r 

for some integer k. In that case, the semiclassical energy E^ s ^ — H(Ik) is a good approximation of the true eigenenergy 
Ek, and the associated semiclassical eigenfunction ip( qs ) = V'Jf ^ fulfills 

U^ s) = exp (-iEi qs) r/h) + 0(h 2 ) . (10) 

It is, at least in principle, possible to improve the above approximation to an arbitrary order in h. Nevertheless, it 
should be born in mind that Eq. (jl(]) docs not necessarily imply that ipk ^ s an approximation of the true eigenfunction 
ipk of U (or H). This becomes particularly relevant for systems that are invariant under some discrete symmetry — 
say, e.g. the inversion P : q i— > —q — which is such that the invariant curve obeying the quantization condition 
Eq. (^]) and its symmetric partner r' fe = PTk are distinct. In such circumstances, the semiclassical wavefunctions 
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ip k qs ^ constructed on T' k will be the symmetric equivalents of ip^^ > an d the corresponding semiclassical energies E^^ 

and E k qs ^ will be exactly degenerate. 

Since P admits only representations of dimension one, there is, however, a priori no reason that the two exact 
eigenergies are degenerate. Classically forbidden processes, that we generically refer to as tunneling events even when 
no potential barriers are explicitly involved, will generally give rise to an exponentially small (in K) coupling matrix 

element a — \ (ip^^ \H \ip k 9S ' > ) |. Using standard WKB methods, this matrix element can be evaluated semiclassically. 
For instance in the case considered above where P is the inversion symmetry relating two invariant curve T/. and T' k , 
one obtains 



T 

where T is the classical period on the torus and 

a = Im 



H 

exp 



K) 



pdq 



r fc -»r; 



(12) 



is the imaginary part of the action integral taken on a path joining Tfe and rj. on their analytical continuation in the 
complex phase space (see in this context also |j35| ). This is illustrated in Fig. |j where we plot the analytic continuation 
of an invariant torus and its symmetric counterpart in the Harper model H = cosp + cosg. 

The projection of H on the subspace generated by ipi an d V'J 9 ^ then reads (with the proper choice of their 
phases) 




Therefore, although the eigenenergies are only slightly shifted with respect to E k qs \ yielding a splitting 5Ek = 2a, 

(and thus an eigenphase splitting 8(pk — 2ra/K), the true eigenstates are not ip^^ and ip k qs ^ but their symmetric 
and antisymmetric linear combinations. Arnold p3| has suggested to call the semiclassical wavefunctions (Eq. (@)) 
quasi-modes to stress that, although they may fulfill the Schrodinger's equation up to an arbitrary order in fi, they 
are not necessarily an approximation of the true eigenstates. Intuitively, this can be seen from the propagation of a 
wavefunction that is initially prepared on one of the tori T^. Although Eq. ([n]) is fulfilled for a single iteration of U , 
the population of the wavefunction will, after a very long time (or a large number of iterations), be fully encountered 

on the symmetric torus T' k , and oscillates between ip k qs ^ and ip k qs ' > with an exponentially long period 2irh/5Ek. 

In the quasi-integrable regime we consider in the following, quasimodes can again be defined, and one can still 
observe tunneling between symmetry related quasimodes which are degenerate at the EBK approximation. We shall 
see, however, that the way the tunneling mechanism takes place is sensibly more complicated than the two-level 
process sketched above in the integrable case. 



B. From integrability to quasi-integrability 



We consider from now on a system with a Hamiltonian which depends on a small parameter e in such a way that 
the dynamics is integrable for e = and non-integrable otherwise. For sufficiently small but finite values of the 
perturbation, the system will display a quasi-integrable dynamics, which more or less means that the classical motion 
is visibly not distinguishable from an integrable one. As stated by the Kolmogorov-Arnold-Moser (K.A.M.) theorem 
(cf. |33|), the phase space of such a near-integrable system is still characterized by dense layers of invariant tori - 
so-called K.A.M. tori - which are slightly deformed with respect to the integrable limit. 

This modification of the phase space structure can be explicitly reconstructed by means of classical perturbation 
theory. Using for instance the Lie transformation method p9| , a (time dependent) canonical transformation (p,q) i— > 
(p(p, q, t), q(p, q, t)) of the phase space variables (p, q) can be defined in such a way that the Hamiltonian is effectively 
time-independent in these new coordinates. This procedure is described in detail in appendix ^ for the special case 
of rapidly driven systems (where e is given by the period of the driving) . Generally, it yields the new Hamiltonian as 
a power series in the perturbation parameter e, which in practice is iteratively calculated up to some maximum order 
n: 



G 



H 6 ~ (p, §) = £ if, (p, g)e ; . (14) 
z=o 

The convergence of this series is in general of asymptotic nature, which means that for any finite e the development 
converges up to some optimal order no and starts diverging beyond. 

As is well known and as was first emphasized by the Poincare Birkhoff theorem (cf. 33 1), the development 

diverges particularly fast in the vicinity of nonlinear resonances. If the frequency of the oscillation generated by He 
- given by fi(J) = dH^/dl in the action angle variables (1,9) of He - is a rational multiple of the frequency 
uj = 2n /t that characterizes the time-periodic perturbation, then even a small strength of the perturbation causes a 
substantial modification of the phase space structure. Except for a stable and an unstable periodic orbit, the resonant 
torus and the tori in its immediate vicinity are broken. At their place, a new regular substructure is appearing which 
is winding around the stable orbit and which manifests within the Poincare surface of section in form of a chain of 
eye-like structures, so-called "resonance islands" (we use this terminology in analogy to mixed regular-chaotic systems 
where they may appear as "islands" of regular motion embedded into a "sea" of chaotic dynamics). This island chain 
is separated from the remaining set of the unbroken K.A.M. tori by a tiny chaotic layer which originates from the 
separatrix structure associated with the unstable fixed point. Compared to the size of the resonances, the extension 
of such chaotic layers is practically negligible if the perturbation is rather small and if overlaps of different resonances 
do not play a role [B6l . 

As a typical example, Fig. |^ shows the phase space portrait of the kicked Harper map in the near-integrable regime 
(t = 1). In comparison with Fig. [|, we see that the phase space structure does not substantially differ from the 
corresponding integrable limit. The most significant modification is in fact the appearance of island chains which are 
induced by nonlinear resonances between the kick periodicity and the free oscillation. 

However, despite the overall regularity of the phase space at that strength of the perturbation, the tunneling process 
is already substantially modified with respect to the integrable case. This was already discussed in the Introduction. 
It is illustrated in Fig. |l| where we show the scaling of a typical tunneling rate with the quantum coarse graining. As 



will be explained in more detail in section IV, we plot here the level splittings (or, more precisely, the difference of 
the evolution operator's eigenphases) between the symmetric and antisymmetric states constructed on the tori shown 
on Fig. | We see that the tunneling rates do not follow the smooth and monotonous decrease with l/H that was 
predicted for integrable systems, but exhibit rather significant fluctuations. Moreover, the tunnel splittings are by 
many orders of magnitude larger than the ones calculated from the integrable approximation ( |l4| ) (dashed line) which 
otherwise reproduces the near-integrable phase space structure quite well. 

These findings are in accordance with the fact that the method of analytic continuation of the phase space tori 
to complex domain, which essentially provided the basis for the semiclassical description of tunneling in integrable 
systems, does not work in the nonintegrable case. It is obvious that the two equivalent tori between which we consider 
tunneling do no longer form a single smooth manifold in complex phase space if the dynamics is not integrable (since 
such a manifold would imply the existence of an additional constant of motion). This alone, however, does not neces- 



sarily disable continuation methods of the kind that was described in section II A. If the manifolds that correspond to 
the analytic continuation of the two equivalent tori happen to intersect under some finite angle somewhere in complex 
phase space, then the respective semiclassical wavefunctions Eq. (^]) can be continued until that intersection line, and 
their splitting can be evaluated by means of their overlap at that line. As has been demonstrated by Wilkinson [B7I , 
this yields essentially the same exponential decrease of the splitting with 1/h as in integrable dynamics, but with a 
different power of h in the algebraic prefactor. 

In reality, however, the analytic continuations of the tori do not meet each other, but are interrupted at their natural 
boundaries, consisting of lines of singularities in complex phase space. This phenomenon has been discussed in detail 
by Greene and Percival J15| for the case of the standard map: by means of the Fourier representation of the K.A.M. 
torus as a function of the angle variable, the location and nature of these singularity lines were analyzed, and it was 
found that the complex tori acquire a fractal-like structure in their vicinity. 

This behavior is qualitatively confirmed for the kicked Harper map. By means of an optimization program which is 
described in Appendix [b|, we are able to identify the complex invariant manifold that corresponds to the continuation 
of a given K.A.M. torus. A typical example of such a manifold is shown in Fig. ^. Although the dynamics is rather 
close to integrability, the K.A.M. torus cannot be continued far away into imaginary domain. In fact, the projection 
of the complex torus to real phase space is restricted to regions far inside the regular K.A.M. tori regime - i.e., far 
away from the chaos border around the separatrix. At this stage of our investigation, we therefore note that the 
concept of analytic continuations does not seem to represent the appropriate framework for the semiclassical study of 
near-integrable tunneling phenomena - which again indicates that in near-integrable systems tunneling proceeds in a 
way that is very different from integrable ones. 
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III. NONLINEAR RESONANCES AND TUNNELING 
A. Effective Hamiltonian in the vicinity of resonances 

As we have seen in Fig. || and Fig. ||, the major modification of phase space that arises when going from integrable 
to near-integrable dynamics is the appearance of substructures due to nonlinear resonances. It is therefore natural 
to ask what would be the influence of these substructures on the tunneling process. In fact, the effect of resonances 
on semiclassically quantized energy levels and eigenstates in a quasi-integrable system has already been discussed 
under a variety of aspects, mostly within the chemical physics literature jl^-|2^]. The approach that we are adopting 
follows more or less the lines of the derivation undertaken by Ozorio de Almeida [^0) : we shall introduce an effective 
integrable Hamiltonian for the dynamics in the vicinity of the resonance and then discuss, in the following section, 
how this Hamiltonian may induce couplings between integrable eigenmodes. 

Most conveniently, the effective integrable Hamiltonian that generates the dynamics in the vicinity of a nonlinear 
resonance is constructed by means of secular perturbation theory [ p9[ . This procedure is illustrated hereafter for the 
particular case of a periodically driven one-degree-of-freedom system. For this purpose, we write the Hamiltonian of 
our weakly perturbed system in the form 

H(I,e,t) = H (I)+e'V(I,6,t) (15) 

where Hq represents the unperturbed Hamiltonian or a suitable integrable approximation of H, obtained e.g. by 
standard classical perturbation theory as a series of the form ( |l4| ) , and (I, 9) are the action-angle variables associated 
with Hq (which implies that H(I,6 + 2n,t) = H(I,6,i)). The perturbative term e'V contains then "all the rest" of 
the Hamiltonian H and is simply given by H (I, 9, t) — Hq(I) for a particular choice of Hq. 

A resonance condition arises whenever the frequency to — 2tt/t of the external driving equals a rational multiple of 
the internal oscillation frequency 0(1) = dH /dI of the system - i.e., 

rQ r:s = slu (16) 

where r, s are coprime positive integers and fl r:s = fl(I r:s ) is the oscillation frequency at the action I r:s at resonance. 
In the vicinity of such a r:s resonance, standard classical perturbation theory diverges rather quickly due to small 
denominators. To avoid this problem, it is convenient to perform a canonical transformation to the frame that 
co-rotates with the angle variable 9 on the resonance. This is done by introducing the new angle variable 

9 = 9 - fl r:s t (17) 

which remains constant, under the time evolution generated by Hq, on the r:s resonance, and varies slowly in its 
vicinity. After the corresponding transformation H i— > fj( r:s ) = H — f2 r:s I of the Hamiltonian (which is necessary 
since the transformation (jlT]) depends explicitly on time), we obtain 

H^ r:a \i,e,t) = H {i)-n r .. s i + e'v {r:s) {i,e,t) (is) 

as new Hamiltonian that describes the time evolution of the new phase space variables (/, 0), with the perturbation 
term 

K< ri '>(j,e,t) = v (i,e + n r:s t,t) . (19) 

Since varies on a time scale that is rather long compared to the periodicity 2tt/uj of the external driving, we can 
now apply adiabatic perturbation theory to the Hamiltonian |29f| and eliminate the explicit time dependence by a 
canonical transformation (I, 0) i— > (1,0) to new, slightly shifted phase space variables (7,0), which is accompanied 
by the transformation i— > fj( r:s ) of the Hamiltonian. In lowest order in the perturbation, this amounts to 
replacing y( r:a ) by its time average over r driving periods (note that V^ r:s \ as defined in Eq. (p"9|), is litr juj periodic 
in t). We thereby obtain the effective time-independent Hamiltonian 



with 



H (r:s) (1, 6) = Ht s) (I) + V {r:s) (I, 6) , 

Ht s \i) ■■= h (i) - n r J , 

V {r:s \l,Q) := — f T e'V {r - s) (I,Q,t)dt. 
rr J 



(20) 

(21) 
(22) 



The new phase space variables (I, 0) are given by 

/ = /+^(/,e,t), ( 23 ) 

e = e + ^(i,e,t) (24) 

where G is, in first order in the perturbation, evaluated as 

G(I,Q,t) = - [ (e'y (r:s) (^e,i') ~V^ s) (lQ))dt'. (25) 
Jo 

Further insight into the properties of the effective Hamiltonian Eq. ( p0| ) is obtained by the Fourier series ansatz 

oo 

e'V(I,e,t)= V Kl (I)exp(i(ke + lujt)) (26) 

kj — — oo 

for the perturbation term in Eq. (jl^) (with V-k,-i — V£{). This yields according to Eq. ( ^2|) 

oo 

VW(I, 0) = V (I) + V r . m (I) cos(mr0 + <p m ) , (27) 



m— 1 



where the real expansion coefficients T^. m and their associated phases (p m are introduced via 

V mr ,- ms = 2 V;. m e^ m . (28) 



We note that, in lowest order in the perturbation, the effective Hamiltonian Eq. ( pp[ ) corresponds to a 27r/r periodic 
function in 0. 

In general, the magnitude of the Fourier coefficients V r . m decreases rather rapidly with m. More precisely, assuming 
the perturbation V to be an analytic function, the decrease of \Vk,i\ with k would be exponential, i.e. 

lim MW^lj = d . (29 ) 



with the exponent d gove rned by the location of the singularities of V. We shall come back in section IIIC, and in 
more detail in section [TV] in the particular case of the kicked Harper model, to the descriptions of these singularities. 
Expanding Hq up to second order around the action I r:s of the r:s resonance, we then obtain 
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= H (Ir:s) + V Q (I r:s ) + (/ 2 /r:s) + Vr.m(Ir:s) COS (mr0 + <p n ) , (30) 

" m=l 

as lowest order expression for the integrable Hamiltonian, with fi = (d 2 Hq / dl 2 )^ 1 (I r:s ) the inverse change of frequency 
with action at the resonance (note that the linear dependence on (7 — I r:s ) is canceled by the £l r:s I term in (pl|)). As 
the first to = 1 term dominates the Fourier expansion of the dependent part, the dynamics generated by Eq. ( |30| ) 
is nearly equivalent to the one of a generalized mathematical pendulum, containing r regions of bound motion within 
< < 2-7T instead of one. At large deviation from resonance ((J — Ir-.s) 2 / (2/i) ^> maxe V 1 ^ 7 "'^ ~ V r .i), the system 
performs a free rotation in that is only marginally perturbed by the presence of t/( r:s ), while in the immediate 
vicinity of (I = I r:s ), librational motion around the (co-rotating) angles Q = 2ln/r (I 6 Z) is also possible. 

This is illustrated in Fig. ^ where we plot the phase space of the kicked Harper map at r = 1 as a function of the 
action-angle variables (J, 9) associated with an integrable approximation of type Eq. (|l4j) . Clearly, we see that the 
structure of the 10:1 resonance chain strongly resembles the one of a pendulum with 10 islands. 



B. Transitions in the generalized pendulum 



The most straightforward way now to quantitatively analyze the effect of the resonance onto the unperturbed 
eigenstate is given by the framework of quantum perturbation theory, directly applied to the effective pendulum 
Hamiltonian (pQ). This requires that quantities like energies, matrix elements, transition rates etc. remain invariant 







under the succession of canonical transformations that leads from (p,q) to (I, 9), which is generally fulfilled in the 
semiclassical regime. We therefore consider now the quantum Hamiltonian 

H = H + V . (31) 

The unperturbed part 

{I -Irs) 2 , N 

H = (32) 



is more or less equivalent to the initial integrable Hamiltonian Hq(I) (see ([15])) near / = I r:a with the term J7 r:s J 
being substracted, and the perturbation which contains the effect of the resonance is given by 

oc 

V = ^ v r.rn cos (mr® + tp m ^J (33) 

m— 1 

(constant terms are omitted throughout). / and being canonically conjugate variables, the action operator is defined 
by 

~ h d 

7 = 7M< (34) 

(with anti periodic boundary conditions in 9 to account for the Maslov indices in the original variable pOJ). The 
unperturbed modes \ipk)i which correspond to the eigenfunctions of Hq, are then given by plane waves in the angle 9 

(e$ k ) = -L cxp (i(k + 1/2)9) , (35) 

VZTT 



with the quantized actions Ik = h(k + 1/2). Their associated energies (with respect to H ) read 

Ek = J-(fc - « r:s ) 2 (36) 

with K r:s + 1/2 = I r:s /h. 

The matrix elements of the perturbation operator V Eq. ( p3| ) within the unperturbed basis are evaluated as 

~ ^ ~ i r 27r , - 00 v 

= ± / V(Oy( k ~ k > e d9 = J] I^5| fc _ fc1 , ror . (37) 

^° m=l 
Hence, within a perturbative approach, the modification of the eigenmode \tpk) reads 

\lpk) = \lpk) + ^2 -Ak,mHk ) (38) 
m^O 

where in first order approximation, the transition amplitudes are given by 



,(1) _ ^|r.m|/2 
Ek — Ek+r 



4 J. = ;'7 ■ (39) 



Second and higher order corrections contain sums over products of type (K-.m/2)/(-Efc — Ek+mr) (V r , m ' /2)/(Ek — 
Ek+( m +m')r) I ^fe+(m+m')r) ■ ^ s a consec i uence : a r:s resonance couples, as expressed by Eq. (p8|), only those unper- 
turbed modes to the state \ipk) the quantum numbers of which differ from k by integer multiples of r. 

The perturbative expansion converges rather fast as long as \V r .i\ <§C |-Eft — Ek± r \ - that is, with Qk = (Ik — E-.s)/ 

rh\Cl k \ > , (40) 

which is well fulfilled as long as the action range spanned by the librational islands is small in front of rh. Due to 
the exponential scaling of V r . m with m (see Eq. (p9|)), the resulting overlap matrix elements \(ipk+rm\4'k}\ decrease 
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in general rapidly with to. Significant admixtures, however, are induced from states \ipk') the quantized actions Ik' 
of which are located on the other side with respect to the pendulum center at I r:s and lie close to the symmetric 
equivalent of Ik- In this case, (ft + k')/2 ~ K r:s , or equivalently, k' = k + rm with m^O such that 

rm ~ 2(k,. :s - ft) , (41) 

and the energies (|36] ) of the states and |^/>fc') become near-degenerate, which strongly enhances their coupling 
with respect to the neighbors \ipk+(m±i)r) ■ Though relatively weak as compared for instance with the \ipk} <-> IV'feir) 
couplings, these transitions across the island chain play a crucial role in the tunneling process. 

This makes it necessary, however, to consider somewhat further the perturbation expansion. Indeed, the exponential 
behavior Eq. (|2^) of the coefficient V r . m makes it a priori not obvious to decide whether, in the evaluation of Ak, m — 
{ipk+mr\^k) , the first order contribution in V r . m dominates the m'th order contribution in V r ,i, since this latter is 
proportional to (V r .i) m , and therefore both terms have an exponential part exp(— md). As we shall see, it turns out 
that the amplitudes (V'fe+mrlV'fc) are dominated by the first order term J^ k m in the limit of small perturbations (at 
fixed h), while for more strongly perturbed systems (or deeper in the semiclassical regime at fixed strength of the 
perturbation) higher order coupling terms may become dominant. 

For this purpose, it is useful to consider in more detail the special case of the exact pendulum dynamics 

(I -I ) 2 

H = - — + V rl cosr6, (42) 

with V r . m — for to > 1. In this case, the coupling from \ipk) to \ipk+rm) is described by perturbation theory of order 
m, which can be straightforwardly evaluated due to the tight-binding structure of the Hamiltonian matrix. As shown 
in appendix one has for fc < K r:s 

ti k+mr m^(§j-) m (43) 

\ h z r z J to!1 (to + Ok) 
in the limit of large to ~ 2(n r:s — ft). Here we introduce 

5 k = 2 Kr:s ~ - m = s k /(hn) - m , (44) 
r 

the equivalent of the energy denominator in terms of quantum numbers, where 

Sk =2 7 rh'^— - (45) 
r 

represents the phase space area that is enclosed between the quantized torus Ik and the center I r[s of the pendulum 
within the angle range < < 2n/r. 

From the semiclassical point of view, the transition from \ipk) to its counterpart on the other side with respect 
to the pendulum center corresponds to a dynamical tunneling process. Unless n r:s is integer or half-integer, this 
tunneling process is, as in the case of a non-symmetric double well, a non-resonant one, which means that the states 
that are connected by tunneling are not quasi-degenerate, but well separated in energy - or, alternatively formulated, 
that quantized tori on one side of the barrier are connected to non-quantized ones on the other side. Under such 
circumstances, only a tiny fraction of the population may be encountered on the forbidden side of the barrier. 

Based on this point of view, we can derive, by means of WKB theory [ ^2|j3^ ], a semiclassical expression for the 
wavefunction within the generalized pendulum, which includes the tunneling component beyond the pendulum center. 



This construction is shown in detail in appendix JCl in the case of the exact pendulum Eq. (42), and can be generalized 
straightforwardly in the more general case Eq. (133) . It yields 



^ C (e) = ' 



'Tk\Q k 

as semiclassical eigenfunction of the state \ipk)i with 



k — fc+mr 
k'~2K r .. s ~k 



(46) 



5 fc (6) = / J k (&)d& (47) 



o 
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the action integral along the torus and 



5e kk , = (2k t , s -(k + k')) = Ek - Ek ' . (48) 

Here, / = 1^(6) parametrizes the quantized torus associated with the excitation k (which naturally implies Sk(2n) = 

2nh(k + 1/2)), Ofc = G>fe(9) symbolizes the time-derivative of the angle variable along the quantized torus, and 
Tfc denotes its period, i.e., the classical propagation time that elapses between = and = 2n. The coupling 
amplitude At is given by 

A T = , (49) 

2sin(sfc/ft) 

where Uk denotes the imaginary part of the action along the complex classical manifold that connects the quantized 
torus with its symmetric counterpart, and Sk is given by ([45|). Interestingly, these two actions cf^, s k and their relation 
to H fully determine the transition rate across the resonance in the semiclassical limit. 

The semiclassical expression ( [46|) is explicitly derived in appendix ^ for the special case of an exact pendulum 
dynamics (|4|). As it is based on the topology of the phase space structure rather than on the explicit form of the 
potential, we expect its validity also in the presence of nonvanishing (but comparatively weak) higher harmonics. The 



case (42) is nevertheless instructive, as it permits an analytic evaluation of the parameters that enter into (ktq). If 




V(0) = V r .i cosr6, we have 

ah = v^vr.i j y —^ri - dx - ( 5 °) 

Assuming V r .\ <C Ek, i.e. that the quantizing torus k is far away from the librational islands of the resonance, one 
can use that /(£>) — Ik + (K-.i/^fc) cos(r0). If furthermore the perturbative condition Eq. (^) applies, we have 
exp(iSk{Q)/fi) — exp(i/c0)(l — iV r .\ sin(r0) / 'rKlk) and 

' exp ( ~S k (e)) ~ &(0) - -^-(^ fc+r (6) - &_ P (0)) . (51) 
\n J 2rnilk 



n\e k \ 

In the regime n r - s — k^>r, this readily gives the first order perturbation Eq. ( |39| ) (with only V r .\ non-zero). Moreover, 
we verify in appendix ^ that the insertion of ( |5l| ) into the semiclassical expression (^6|) of the k' component recovers 
the quantum transition amplitude d43| ) in the limit |<5fc| <C 1 (or, more precisely, in the limit \Vr.i I -C \Ek ~ Ek>\ <C 
\Ek — Ek'± r \), where only one quantum state from the other side is significantly coupled |p9[ . The semiclassical 
expression becomes particularly useful when the condition ( |4C| ) for quantum perturbation theory does not hold any 
more. 



C. Determination of the coupling strength 

The description of the local dynamics near a r:s resonance by the Hamiltonian Eq. ( |30| ) gives rise to a mechanism by 
which the quasimodes located on opposite sides of the resonance are coupled. This will constitute the basic ingredient 
to the global tunneling mechanism which we shall develop in the next subsection. To allow for a quantitative prediction 
of the associated transition rates, it is necessary, however, to specify how the parameters V r . m and ip m that enter 
into the expression of H r:s can be computed in practice. The purpose of this subsection is to show how this can be 
done from the classical motion near the resonance. We shall furthermore discuss some qualitative properties of these 
quantities, in particular the asymptotic behavior of the V r . m for large m. 

The only slight technical difficulty we shall need to address here is due to the fact that we consider maps. More 
precisely, the integrable Hamiltonian H^ n \p, q) Eq. (|l4|) has been introduced in such a way that the map T it generates 
is the same as T up to 0(e n ) corrections. In other words g^'™' — , where is the hamiltonian flow generated by 
Ti. = H or H^ n \ However, nothing imposes a priori that for intermediate times < t < t, g^ { '(g,p) ~ gp(q,p) up 
to order n corrections. As a consequence, the original Hamiltonian H = H(I, 8, t) is not well approximated by the 
time-independent expression H^ n \l). 

Starting from the action-angle coordinate X = (/, 9) of H^ n \p, g), we shall therefore first need to define a period- 
ically time-dependent coordinate system X b {X 1 t) = (I b ,9 b ) such that in these new coordinates, the kicked Harper 
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Hamiltonian is well approximated by £f(")(J rh ) for all times, up to smal l corrections that we can then deal with by 
using the standard secular perturbation theory described in section [II A| . We are thus looking for a periodically time 
dependent canonical coordinates transformation 7 t b : X \— > ~f b (X) — X b (X, t) such that 

X b (X,r) = X b {X,Q) = X , 
and H b (I b ,0 b ,t) = H(l b ) + O{e no ) . 

A way to fulfill these constraints is to define ^ b as 



It 



= g?og H t , for 0<<<r, (52) 



for < t < t, and by periodicity for the rest of the real time axis, where g^ symbolizes the Hamiltonian flow over 
time t generated by the Hamiltonian Ti. = H or H^ n \ The following scheme 

X gf(X) 

7o b I T (7?)" 1 (53) 

X b (X,0)=X gf(X b {X,0)) 

H 

illustrates why the motion under the Hamiltonian H in the original variables X is equivalent to the one generated by 
H<ri in the X b variables, for t G [0,t[. 

The transformation j b is, by explicit construction, periodic in time. However, it is in general not continuous at 
t = nr, as a consequence of the fact that T does not perfectly approximate T . The complete definition of the new 
Hamiltonian H b requires therefore to introduce a perturbation term e'V b (I b , b , t) which becomes active only at t = nr 
and which accomplishes the final "jump" from X b (X, t — * r) to X b (X, r). One therefore has 

H b (I b , 9 b , t) = H(I b ) + e'V b {I b , 6 b , t), (54) 

with 

00 

e'V b (I b ) e b ,t)= S^(t-nr)R(I b ,6 b ,t). (55) 



e' is the strength of the perturbation induced by V b (I b ,9 b ,t) and corresponds to the accuracy of the integrable 
approximation of H. From a strictly formal point of view this strength is of order 0(e n ). This scaling, however, 
applies only to contributions that are analytic in e (e.g., a global deformation of the K.A.M. tori) and does not take 
into account non-analytical contributions (e.g., of the form C\ exp {—C^jt)) which result from the vicinity of nonlinear 
resonances. 8^ is a Dirac distribution that, for consistency, we need to consider as being smeared on the interval 
[t - ,t], with r~ — > r. (In practice, we take 5^'(t) = l/(r— t~) if t S [t'^t], an zero elsewhere.) R(I b ,O b ,t) is a 
time-periodic function with period r. 

A natural interpretation of what R(I b ,8 b ,t) is can be obtained by integrating Hamilton's equations of motion 
associated with H b (X b ,t) from to r. This yields 

g?(X,T)-g?(X,r) = f JVR(X b (t), t) dt (56) 

T-T J T - 

where X = ( and X b (t) is the path that relates g^(X,r) to g^(X,r). Notice that these equations would be 
inconsistent without a time dependence for R. However, as the path X b (t) is of typical size e', one can rewrite 
perturbatively (^) as 

I?(X) - JfpQ = -&-(I?(X),9?(X)) + 0{{e'f) , (57) 

0?(X)-0?(X) = ?^L{I?{X),6?{X)) + 0{{e'Y) , (58) 

where (R) denotes the time average of R between t~ and r. We recognize that (R) is, in first order in e', the generating 
function of the canonical transformation 
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g?(X,T)^g?(X,r) (59) 
that is, of the difference between the map and the motion of its integrable approximation during a time r. 



Within the X b variables, we can now apply the standard secular perturbation theory described in Sec. IIIAj. We 
obtain in this way 

H r:s (I, 6, t) = Ht s) (I) + V^ S \I, 9) + 0(e' 2 ). (60) 

with 

r-1 



V r . m e lVm = — I exp(-irme){R)(I r:B ,0)d9 . (63) 



V r . m e tipm = — — / exp(-irm9)6I r:s (9)d9 . (64) 



9) = — V(i?)(I, + Q r:s n T ) . (61) 
rr * — ' 

71 = 

The Fourier coefficients of the averaged perturbing potential 

oo 

(I, Q) = J2 V r.m COs(rmO + ip m ) (62) 

m=0 

can, with 9 ~ 9 and I ~ I c_ 7, :s , then be written as 

<-27T 

This transforms after integration by parts into 

r-1-K 

mmrr 

Here, 5I r:s (8) is defined by 

5I r:s (0) = I^(I r:s ,9)~I r .. s (65) 

where /( _1 )(7, 6*) symbolizes the action variable that is obtained by applying the inverse Poincare map T _1 to (J, 0) 
(or alternatively, the backward propagation with H from time t = r to t = 0). Eq. (|64|) therefore provides a 
convenient way to obtain the numerical value of the coefficients V r . m , which is based only on the propagation of 
classical trajectories. 

The effect of averaging out the time dependence on the integrable contributions of V leads to the ^-independent 
coefficient Vo which is of order e n . On the other hand, the other coefficients V r . m with m > correspond to the 
non-integrable effect of the resonances, and therefore their magnitude is not simply proportional to e" (we should 
actually expect them to be essentially independent of n, in some range near the optimal value rio). As the V r . m result 
from the Fourier integrals of SI r[s (9), their scaling with m can be inferred from the analytical structure of 6I r:s (9). 
Assuming 5I r:s (9) to be an analytic function in 9, the line of integration in Eq. j6^ ) can be displaced into the negative 
imaginary direction of 9 (for m > 0), where it g ives a vanishing contribution due to the exponentially small prefactor. 
As a consequence, the Fourier integral Eq. (J6J) is entirely described by the singularities 9^ — fi r:s t^ of SI r - s (8) in the 
complex 9 domain, and will, for large m, be dominated by the contribution of the singularity that is closest to the 
real axis (see in this context also p3f ). 

The calculation of SI r :s(9) involves in practice three steps. The first one is to determine the coordinate (p,q) 
of the point (I r - S ,9) under consideration. The second one is to apply the map T _1 to (p,q), and the last one is 
to determine the action coordinate I of the resulting point. In general, these two latter steps should not involve 
any singularity: the map T, the function H^ n \p,q), as well as the function 1(H) will usually be analytical. As a 
consequence, the singularities of 8I r - s (9) should be the one of the torus T r:s , that is the complex angles 9 such that 
(p, q) = (j>(I r:s ,9),q(I r:s ,9)) lies at infinity. This corresponds to trajectories which, starting from (I r:s ,9 = 0) on the 
real torus, go to infinity in a finite complex time under the dynamics of H^ n ' . 

One can therefore write in the asymptotic regime 

V r .m « K% e -( !) " i ™)~,m> 0, (66) 

with tf° s the imaginary part of the time to reach the closest singularity, and where V^ and 7 characterize the behavior 
of (R) near the singularity. (If (R) was a meromorphic function, 7 would be the degree of the pole, and V^ the 
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corresponding residue.) We would like to stress here is that there are two sources of smallness in this expression. 
One is the exponential dependence in to, which is entirely controlled by the dynamics of the integrable approximation 
(t^° is determined by H^). In the semiclassical limit, this will give rise to an exponential dependence in h, since 
one should use m = AI/rh to connect two tori differing by an action AI. The other parameters characterizing the 
asymptotic behavior of the V r , m , namely 7 and V^ , depend on the complete dynamics of the perturbed system and 
contain in particular the perturbation parameter e. To have a well defined classical perturbation expansion, and in 
particular for the first order secular perturbation approximation we have used to be valid, the corresponding (rm)^V^ 
terms should be small on the classical scale, although not exponentially. We shall always assume the perturbation 
parameter e to be small enough for this property to hold. 

In addition, The general scaling behavior Eq. ( |66| ) has consequences for the quantum perturbative expansion to 

evaluate the overlap (ijjk+rm\'4>k}, and determines up to which order this expansion should be done. To illustrate 
this, let us consider for a particular r:s resonance the coupling between two states that are symmetrically located on 



opposite sides with respect to the resonance (i.e. such that Eq. ( |41[ ) holds). The second order correction to Eq. (|39| ) 
reads 

l7 I / \ 'r.m \ ^ K-.miK-.mj lai\ 

(V>fc+™#fc> = ijrp p r + 2^ JTp W uf W T ■ (67 ) 

If Eq. ( |66|) applies, we see that the condition for the second order term to be smaller than the first order one does 
not involve the exponential, but that for each toi in the sum {m,irn2r/m) 1 V^ a /{Ek — Ek+ mi r) <C 1. For a given 
value of the perturbation parameter e and at fixed h, such a condition may very well be fulfilled. However, in the 
semiclassical limit h — > with fixed to - and in practice to will always be more or less fixed when the system undergoes 
the transition over a particular r:s resonance (as will be discussed in the following subsection) - the denominator 
(Ek — Ek+ mi r) goes to zero (being bounded by Ek — E r:s ~ h 2 (m 2 r 2 /8/i)). Therefore, as one goes deeper in the 
semiclassical regime, the second order term will eventually dominate over the first one. 

In the same way, one can see that assuming Eq. (|6^), the condition for the m th order term Eq. (fl3| ) to be larger 
than the first order one is that 

» 1 • (68) 

In the semiclassical limit (with fixed to), this condition will eventually be reached at one point. 

As a consequence, we see that, assuming V^ to be small on a classical scale, a first order quantum perturbative 
treatment will be valid for moderately small values of h, but higher order should be taken into account as h — > 0. 
Note that this is not incompatible with the fact that the quantum perturbation development is convergent, since the 
condition Eq. (|40| ) for its validity involves the exponential term exp(— rn r:s t^) which can be extremely small, especially 
for high-order resonances with r> 1. Very far in the semiclassical regime (or for small r), quantum perturbation 
theory might nevertheless fail at some point, in which case the it would become necessary to resort to semiclassical 
expressions such as Eq. (ff6|). 

We see that considering the analytical structure of the function SI r:s (9), is important to decide what term in 
the perturbation expansion will be the dominating one, as well as, as we will see in the next section, what is the 
dominating mechanism in the tunneling 'process. This should be reconciled with the fact that the analytical structure 
of the invariant tori may sensitively depend on the precise choice of H^ n \ and in particular on the degree n of the 



integrable approximation. We shall come back on this issue in section IV 



D. Mechanism of resonance assisted tunneling 

In the previous sections, we have examined in detail the characteristics of couplings that are locally induced by 
the presence of a nonlinear resonance. We shall now see how these couplings can be combined at a larger scale to 
form a global mechanism of tunneling for quasi-integrable systems. Furthermore, we analyze why, and under which 
condition, this mechanism is th e dominating one. 



As we have seen in Sec. Ill A, the dynamics near a r:s resonance is locally described by a Hamiltonian of the form 

(T — T ^ 00 

H ir-,) = I + V™ cos (mr9 + <p m ) , (69) 

^ m—l 

where the parameters and tp m can be computed with Eq. ( |6^) through the propagation of classical trajectories. 
Furthermore, when discussing the order of magnitude of the various terms, we shall assume that the asymptotic 
expression derived in the last section can be used, and thus that 
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Km ~ (mryV£ s exp(-JW&mr) , (70) 

where ft r - s = suijr is the angular frequency of the integrable torus at the r:s resonance, tf° s > is the imaginary part 
of the classical time to reach the closest singularity of the analytic continuation of r r:s into complex phase space, and 
7 and V^ :s characterize the generating function (R) near the singularity. 

Let us, to start with, consider the unpert urbe d Hamiltonian H(e = 0), and one of its quasi-modes tp k qs ^ built on 



an invariant torus IV As discussed in Sec. II A, the symmetry P of our system is assumed such that Yk exhibits a 



symmetric, but distinct equivalent Y' k on which one can build another quasi-mode ipu qs ^ analogous to ip^u^ ■ and 
r' fc are separated in phase space by a separatrix S. The true eigenstates of the evolution operator U correspond to the 

symmetric and anti-symmetric linear combination of ipi^ and ipS qs ^ , the eigenphases of which differ by the splitting 
<5</>fc- The semiclassical expression of the splitting is given by 2ar/h, where a is defined by Eq. (pd[). 

Now we increase e and follow the adiabatic evolution of the eigenmodes of U , which can be considered as pertur- 
bations of the quasi-modes ipk!^ associated with the integrable approximation H(e). At some point, a resonance r:s 
grows significantly as compared to h and couples ip^^ to some V'i+rm' with m > 0. Since the torus Yk+rm is located 
closer to the separatrix than Yk, tp^+lm exhibits a slower exponential decrease in the forbidden domain than ifi k qs ^ '. 

As a consequence, if the strength of the coupling between ip k qs ^ and V'i+rra 1S n °t too small, the admixture of this 
latter component will eventually dominate the behavior of the perturbed quasi-mode in phase space regions close to 
the separatrix, and thereby determines the eigenphase splitting 54>k between the symmetric and the antisymmetric 
linear combinations of the quasi-modes. One obtains in this way a splitting 



^, s J 2 ^exp(-^), (71) 

J-k+rm » ' 



where Tk+ rm is the classical period of the torus rfc +rm , ak+rm is the imaginary part of the classical action along a 
complex trajectory relating Yk+rm to its symmetric counterpart, and A k ' s m represents the coupling amplitude ( |3q ) 
associated with the resonance. For sake of clarity, we shall consider below the case where A k s m can be approximated 
by the first order expression Vf^/ '(E k ' s — E k ? rm ) with E k :s — Ek — Vl r . s Ik- Our argumentation, however, does not 
rely on this precise form. 

To compare the relative effectiveness of the above "resonance assisted" mechanism with respect to the direct 
(integrable-like) one, we use Eq. (^) and obtain that 

* a (rm) iF-Bsu exp v ' s ) ■ (72) 

If, for a moment, we just compare the exponential factors of the above expression with the one, exp (— <Jk/fi), of the 
direct tunneling mechanism, we see that the condition for the resonance-assisted one to be dominant would be that 

2Q r ; S tf° s (Ik+rm — Ik) < Cfc — CTk+rm ■ (73) 

Now, in the semiclassical regime, one can assume Ik, Ik+rm and I r:s classically close, and thus T . :s (J/ c + I . m — Ik) — 
AE = Ek+rm — Ek- In the same way, ak+rm — &k — AE ■ (da / 'dE) r - s , with (da / dE) r:s = t^. s the imaginary part 
of the time needed to follow the complex path V from the resonant torus T r:s to its symmetric counterpart Y' r . s , on 
which the action a is computed. As a consequence Eq. ( [73] ) reads 

4% < % , (74) 

or in other words, that the imaginary part of the time needed to reach the closest singularity should be smaller 
than half the imaginary part of the time required to go from one torus to its symmetric partner. This condition is 
necessarily fulfilled, as can be seen from propagating V under for real time t £ [0, t$? s ]. (gf 1 ' V necessarily 

encounters at least one singularity of r,. :s , and by symmetry, one of these singularities necessarily fulfill Eq. ( f74|).) As 



a consequence, the resonance-assisted mechanism will always dominate the "regular" tunneling process (Sec. II A) in 
the semiclassical limit. 

Considering now the prefactor, the energy denominator in Eq. ( |7l| ) will make it favorable to connect i^k*^ to a 

quasi-mode V'i+Lr sucn tnat ^k' s 1S am i° s t (i-e. up to a difference of order Sl r:s /i) degenerate with E k + rm , which 
implies that r r:s should lie at mid distance between the tori Yk and Yk+mr- Note that this is the case not only if 
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the first order approximation of A r k ' s m is used, as in Eq. (|72[), but also if higher order terms of the perturbation are 
included, or if the semiclassical expression Eq. ( pi| ) is used. 

For small ft, Tk+ m r will then a priori not be close to the separatrix. However, nothing prevents from making use of 
couplings via other r':s' resonances in order to gradually approach the vicinity of the separatrix 5. In this way, tp^ s ^ 

can eventually be connected to a quasi-mode ipw^ the action of which is only a few ft smaller than the action Ig 
of the separatrix and from where "regular" tunneling takes place with a rather large rate. 

Using successively the resonances ri'.si, r2'-S2, ■■■ r^-.SN, which we assume to appear in ascending order (i.e., 



< 



< I r 



, ) , the resulting expression for the splitting is then 
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k N 7 m N 
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(75) 



Here, k = k\ , k 2 = k± +toiTi, k$ 



m 2 r 2 , 



fcjy + mN^N denote the quantum numbers of the intermediate 



quasi-modes that are involved in the coupling scheme. The rrii are always chosen such that the denominators E r k 



E. 



ki-\-rim-, 



are quasi-degenerate, which means that the tori and 



should be almost symmetric with respect 



to r rj 

In the particular case where the semiclassical expression Eq. (46) can be used for the amplitudes A™ m , one obtains 
the expression 
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(76) 



with ek i _ 1 k i defined by Eq. (E3) and 
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(77) 



(see Eq. (Eq)) the phase space area that is enclosed between and r ri:Si within the angle range < 8 < 2n/ri 



Here a 



hi 



represents the imaginary action along the complex path that connects, within the effective pendulum 



Hamiltonian Eq. (|69|), the perturbed torus with its counterpart on the other side of the r^:Si resonance. The 
overall picture is illustrated in Fig. |^. Within the corresponding secular perturbation approximation, each resonance 
provides a complex path allowing to join a torus on one side of the resonance to its symmetric counterpart on the 
other side. The global "resonance assisted tunneling" mechanism we propose consists in following this succession 
of complex paths across a series of resonances to reach the neighborhood of the separatrix. From there, "regular" 
tunneling occurs into the symmetric island. 

In the above scenario, it remains to decide what is the most effective chain of resonances to be used in the tunneling 
process. For a given value of ft, one constraint is naturally that a r:s resonance needs to be taken into account 
only if the area 2itls enclosed by the separatrix is larger than 2-Krh. For relatively large ft, this can leave only a 
moderate amount of possibilities. Deep in the semiclassical regime however, a large number of resonance will fulfill 
this condition, and it is useful to design a guiding principle on which one to use. 

For the sake of clarity, we shall address this question again under the hypothesis that the first order ap- 
proximation Eq. ( |39| ) can be used for the transition amplitudes A k s m . With Eq. (JTTj), we can therefore write 
A r k ; S m - B £!m exp(-ft r: 4? s rm), with 



yps 
so 



= (rm)T._ 

^k+rrn 

With this notation, we have, for a given choice of resonances r\:si, ■ ■ •, r^:sN, 



(78) 
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gr N :s N 
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2 2r f 1 . 
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(79) 



with, noting rivrii — A/^ i:Si )/ft, 
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J kk > = ^Vu-.AUA 1 ^ (so) 



h 
»=i 



rt(i)t^(i)di 



The exponential term exp (—(Jkk' + 0V ) / ft) will therefore not depend too much on the precise choice of the sequence 
of resonances used to connect If. and I' k . On the other hand, using a large amount of resonances between and I' k , 
each of them inducing only a small change in the action, will have a tendency to increase the number of B k s m terms 
in the prefactor of 8<f>k- Hence, the choice of the dominant path will in the end depend on the magnitude of these 
coefficients. 

The energy denominator of B k s m is in general of the order of AE r:s — rSA/' r:s '/2/j. Therefore, the condition that 
it is favorable to introduce a r:s resonance into the coupling path reads 

&k,m * h2r S " » 1 (82) 

Even though the parameter V^ is usually small when the system is close to integrability, this condition will eventually 
be met if one goes high enough into the semiclassical regime. This can be interpreted as an upper bound for h above 
which it is impossible to "resolve" the r:s resonance. We stress though that this criterion is not directly related to 
the size of the islands of the resonance. Indeed, this latter quantity is proportional to y/V™ and involves therefore a 
factor exp(— Q, r - S tf.° s r) which can be extremely small for large r. V^ s should, on the other hand, smoothly depend on 



0,/uj — s/r (except for symmetry considerations, see Sec. IV) and might not be very different between, say, Q/uj = 1/2 
and fl/u) = 5/11. Note finally that using higher order terms in the perturbative calculation of A k s m can only make 
the condition Eq. ( |S2| ) valid sooner in the semiclassical regime. 

In any case, the above consideration implies that as long as Eq. ( |82| ) applies for all resonances that are used in 
the coupling sequence, and as long as the final torus is kept as close as possible to the separatrix, a large number of 
small steps will generally be favored with respect to a small number of large ones, since each extra step imply the 
multiplication by a large prefactor (which, roughly speaking, comes from a small energy denominator). This justifies 
a posteriori the use of a local description of the dynamics, since the dominant tunneling mechanism takes into account 
the effect of a resonance only in its close neighborhood. 
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IV. RESONANCE ASSISTED TUNNELING IN THE KICKED HARPER MODEL 



After the general discussion in the previous sections, we shall now consider in more detail a particular system in 
the nearly integrable regime, namely the kicked Harper model. The purpose of this section will be to check, for this 
particular case, the accuracy of the final semiclassical expression Eq. (|75|) - which gives the splitting between the 
quasi-energy of a symmetric Floquet mode and its antisymmetric counterpart - as well as to verify the degree of 
validity of the various hypotheses that were made along the way of its derivation. 

The classical dynamics of the kicked Harper is governed by the hamiltonian Eq. (Q) , yielding the stroboscopic map 
Eq. (||). In the limit r — > 0, this dynamics is equivalent to the one generated by the time-independent (integrable) 
Harper hamiltonian 

H (p,q) = cos p + cos q . (83) 

In Eq. (Q), r is thus both the period of the kick and the perturbation parameter (i.e. e = r). 
Quantum mechanically, the map Eq. (^) can be associated with the evolution operator 

f T ( iTCOsp\ ( ircosq\ 

U = exp exp — — 1 . (84) 



h ) 1 \ h 

The 2n periodicity in q and p makes the quantum treatment of the kicked Harper particularly easy if 

(85) 

with integer TV > 0. For these particular values of K, the eigenfunctions i/j of U can be written as Bloch functions in 
both position and momentum - i.e., 

(q + 2tt|V) = cxp(iip q )(q\ip) , 
(p + 2n\i') = cxp(i(p p )(p\i/j) 

for some pair of Bloch phases < (p q ,ip p < 27r, where \q) and \p) denote the eigenfunctions of the position and 
momentum operator, respectively. For each pair (ip q , tp p ) of Bloch phases, the corresponding subspace of the Hilbert 
space is finite dimensional and contains N linearly independent wave functions, spanned, e.g., by the basis states 

oo 

\n) = ex P(*( z + n/N)cp q ) \Q = ( 27m + <P P )/N + 27r (86) 

/ — — oo 

for < n < N f2~4}] . The eigenvectors \ipk) of U and their eigenphases cf>k can therefore be computed up to numerical 
(quadruple) precision, by diagonalizing the NxN matrix (n\U\n'). In the following, we shall consider only the two 
pairs (0, 0) and (0, it) of Bloch phases, corresponding to periodic boundary conditions in momentum, and periodic or 
anti-periodic boundary conditions in position. This choice is equivalent to restricting p to the interval [— 7r, 7t] and q 
to the interval [— vr, 37r] with periodic boundary conditions (see Fig. ||), and to consider the even and odd symmetry 
classes with respect to the inversion q — > —q. 

The calculation of the integrable approximation H for the kicked Harper is performed straightforwardly by applying 
the formalism of Appendix One obtains for instance as zeroth order coefficient the Harper Hamiltonian Eq. pS] ) , 



and Eq. (A2E) for the approximation of order three (recall that e = t). In principle, one may co nstru ct up to 

orders as high as n — 20 fairly easily with symbolic programs such as MAPLE. As mentioned in Sec. [II B| , however, the 
series Eq. @ of tends to re-diverge beyond an optimal order no, which, for r = 1., is generally found around 
no ~ 6. This is illustrated in Fig. |[ For various orders n of the integrable approximation, 40 randomly distributed 
initial phase space points have been propagated during a given time by means of the kicked Harper map as well as 
by its integrable approximation H^ n \ and the distance in phase space between the two resulting sets of final points 
is plotted as a function of n, yielding a minimum at rather moderate values (n ~ 6 in this particular example). We 
shall therefore mainly use with n = 6 in the following. 



19 



A. Resonances parameters 



Fig. |]and Fig. || compare the phase space portraits of the kicked Harper and of its integrable (6th order) approxima- 
tion in the near-integrable regime at t — 1.0. In fact, one observes that the only significant difference between the two 
Poincare sections is the presence of the resonances. One may further note the relative importance of r:l resonances 
with r = 10 and 14 as compared to the 8:1 and the 12:1 resonances (the absence of resonances with odd r is an 
obvious consequence of the rectangular symmetry of the kicked Harper). As a matter of fact, these latter resonances, 
with r a multiple of 4, are rather weakly developed at r = 1 and systematically exhibit 2r (instead of r) islands in 
the Poincare surface of section. We conjecture that this behavior is a consequence of the initial square symmetry of 
the Harper Hamiltonian, which is still relevant for small values of r. As the period in the center of the regular region 
is already larger than 6 and monotonously increases when moving towards the separatrix, r:s resonances with r < 6 
do not exist at r = 1. 

To obtain a quantitative prediction for the tunneling rates, it is necessary to characterize the resonances through 
the Fourier coefficients V£{£. This is done in practice by a direct application of Eq. (|6J), i.e. by Fourier transforming 
the function SI r:s (9) = 5I rts (p(I r:s , 9), q(I rts , 9)) where I r:s is the action of the resonant torus T r:s . On this torus, the 
angle variable is given by 9 — O r:s £, with Q r . s = dH/dI{I TXS ) = 27rs/(rr). For a given 9, 6I r:s is computed through 
the following successive steps: i) Choose once for all a reference point {po 7 qo) on the resonant torus T r . s of H^ n °>. 
ii) Propagate (po;<?o) under H( n °) dynamics during the time t = [r / s){8 /2-k)t. iii) Apply the time reverse of the 
Poincare map Eq. (|2j) on the resulting point, iv) Compute the difference between the action of this iterated point 
and the action I r:s of T r[s . The values obtained in this way for the 10:1 resonance are plotted on Fig. § for various 
orders n of the integrable approximation, showing that for 3 < n < 6 the coefficients VI £ do not depend sensitively 
on n. Also shown in this figure are the values obtained by the method introduced in p3j, which is based on a Fourier 
analysis of the (pseudo-)separatrix structure that is associated with the resonance. 

Within our setting for the kicked Harper, the tunnel splitting is defined as the difference 

S<t>k = \MVi = 0) - <t>k(<Pi = tt)|. (87) 

As already stated, the exact quantum values of S<f>k can be calculated up to numerical precision. Using the coefficients 
obtained in the above way, as well as the unperturbed energies Ef., the periods Tk and the tunneling actions 
ct/c which are straightforwardly calculated from the integrable approximation H^ n > of the kicked Harper, these exact 
splittings can be compared with the ones derived from our semiclassical expression Eq. ( |75"| ) based on the resonance- 
assisted tunneling mechanism. 

Before perfor ming this comp arison, let us first verify that the qualitative description of the tunneling mechanism 



we gave in Sec. 1IIB and HID actually applies in this particular example. To start with, we can check that all the 
resonance involved in the tunneling process are well within the quantum perturbative regime. Indeed, for the value of 
the perturbation parameter we consider, r = 1, the largest Fourier coefficients for the resonances coming into play are 
Vf 6 :1 ~ 9.0 • 10~ 7 , (as already stated, the 8:1 resonance exhibits 16 islands), V^ 1 ~ 2.5 • 10~ 4 , V$ :1 ~ 9 • 10" 4 , while, 
in the range of h we consider, the energy difference between quasi-degenerate states with respect to the resonance 
is typically of the order of ~ 10 -2 . Furthermore, taking into account the actual values of the V£ s r we observe that 
as h = 2it/N gets smaller, higher orders of the quantum perturbation theory become dominant in the calculation 
of the transition amplitudes Ap." L. This can be specifically verified for the 10:1 resonance: For this resonance, the 
k — > k + 20 transitions are of order one - i.e., are dominated by the first-order perturbative coupling terms - for 
N < 38, but involve perturbation theory of order two for N > 38. Similarly, we find that the k — ► k + 30 transitions 
are of order one for N 5 38, of order two for 38 5 N 5 127 and involve higher terms beyond (k — > k + 10 transitions 
are, of course, always of order one). Effectively, one finds here the (possibly unusual) situation which will generally 
be encountered in the semiclassical limit - namely that the lowest order terms of the perturbative expansion (which 
converges nevertheless well) are not the dominating ones. 

Figures and ^ show for a varying value of ft, i.e., a varying total number N = 2n/h of states, the eigenphase 
splittings of the eigenmode of U that corresponds to a fixed classical torus, with action I = 7r/4 in Fig. [I] and with 
action / = tt/6 in Fig. [Io|. Evidently, these splittings can be calculated only for particular values of N, namely for 
N = 4(2fc + 1) and N — 6(2fc + 1) with k = 0, 1, . . . in Fig. [l] and Fig. [l(], respectively, for which this torus is selected 
by semiclassical quantization and supports the fcth excited quasi-mode. In both cases, the perturbation parameter r 
equals 1.0. The resonance involved are the 8:1, 10:1, and 14:1. We observe that the agreement between the quantum 
and semiclassical results is extremely nice. For the moderately small values of h that we consider, it is possible to 
try all the possible coupling paths that participate at the tunneling process, and in Figs. [I] and [h], the semiclassical 
prediction is obtained by summing up all these contributions. However, as shown on the lower panel of Fig. |l^, where 
the action coordinates of the intermediate states that participate at the dominant tunneling path arc displayed, we 
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see here that, as discussed at the end of section [II D, this dominant path is always such that the number of steps is 
as large as possible, taking into account the contraints due to h. 

Finally we show on Fig. |ll| a comparison, for a fixed value of h and a variable initial torus, between the exact 
quantum mechanical splitting and the one calculated from the expression corresponding to integrable tunneling, with 
no resonance coupling. We observe on this figure that, although the two curves strongly differ in the interior of the 
regular region, they match perfectly as one gets close to the separatrix. This shows that the presence of the separatrix 
does not introduce any additional effect (e.g. from a small chaotic layer) to the tunneling mechanism. 



B. Singularities of the invariant manifold of the integrable approximation 

In addition to the numerical values of the coefficients V££, needed to obtain quantitative prediction for the tunneling 
rates, a qualitative understanding of their behaviour, and in particular their asymptotic properties for large m, is, as 



seen for instance in section HID, also required to guaranty that the tunneling mechanism we propose is indeed the 
dominating one. Since the are proportional to the Fourier coefficients of the function 5I r ; S (9), their asymptotic 
behaviour is related to the singularities of this function, for complex values of the angle 9. 

Let us consider, more generally, for fixed values of the energy E and the order n of the integrable approximation, 
the invariant manifold T of H^ n \ defined by the equation H^(p, q) = E and characterized by the angular frequency 
SIq. Let a function f(9) be defined on T as 

f(O) = f(p(0),q(6)), (88) 

where f(j>, q) is an entire function of the phase space variables. As a consequence, the singularities of f(9) are the ones 
of (p, q)(9). What we therefore need to study are the singularities of the analytic continuation of (p, q)(0) for complex 
angles 9. Due to the linear relation between 9 and t, this analytic continuation is straightforwardly constructed by 
propagation (under H^ n ') of some real initial point (po, go) on I\ taken as the origin of the angle axis, over complex 
time t. A singularity of (p, q)(9) is an angle 9^ such that for the time ft = the point g^(po, qo) goes to infinity. 
Note that because of the existence of these singularities, (p,q)(9) actually depend not only on the final time 9/CIq, 
but also on the homotopy class of the path joining t — to 9/Qq in complex the time plane. In other words, (p, q)(9) 
is a priori a multivalued function of 9. 

To search for the singularities of (p, q)(9), the first step will consist in finding asymptotic expression describing the 
manifold H^(p,q) = E when the imaginary part of p and/or q goes to infinity. For this purpose, we introduce the 
variables 

A = cxp(ip) , (89) 
Y = cxp(ig) . (90) 

In these new variables, the integrable approximation of the kicked Harper Hamiltonian takes the polynomial form 

+(n+i) 

M = -(™+l) 

with known real coefficient a\™ . For for instance, the non-zero coefficients are a^\ = = a < ^ > \ = = 1/2. 

The manifold T is invariant under the symmetries so ■ (A i— > 1/A, Y i— > l/Y), sa : (A «-> Y), and s : (A t— > 
1/A, Y i— > l/Y). Moreover, one can check easily that if H^[t](X,Y) = E, then H^[—t](1/X, Y) = E. We shall 
call s T this transformation, although this is not properly speaking a symmetry of T. The asymptotic regions of T - 
i.e., the neighborhood of points at infinity on T - can be obtained by application of one of the above transformation 
from one region such that Im[p] — > +oo, (i.e. A — > ) and Im[g] is either bounded or goes to +00 (i.e. Y bounded). 
For such regions, one can assume an asymptotic expression of the form 

y«(A) = 7o « + 7l «A + 7 | A 2 + • • • . (92) 

where £ label the asymptotic region. Introducing Eq. (|9^) in the expression Eq. ( |9l| ) of the Hamiltonian to solve the 
equation H^(X,Y) = E yields a series of polynomial equations for the coefficients 7, , which can be solved order by 
order to determine successively 7(5, 7^, 7I, etc.. Again for the zeroth order Hamiltonian H^°\ the set of equations 
obtained in this way are 
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7o = 

2 (7i + 1 + To) = E la 
7 2 /2 + 71 70 + 70/2 = £71 

yielding 70 = 0, 71 = — 1, 72 = — 2E, ■ ■ ■. In other words, for small X, the manifold defined by the implicit expression 
H(°'(X,Y) = E admits the explicit asymptotic expression 

Y(X) = -X - 2EX 2 + ■■■ . (93) 

Using the above equation with X small enough allows to find a point with a large imaginary part for p, such 
that H w (p,q) is very close to E. This point can be brought back to the energy E by following the gradient of the 
Hamiltonian, giving a point (jfi, q^) on the H^°'(p, q) — E manifold and in the asymptotic region of large Im[p]. From 
this point, we integrate Hamilton's equations of motion choosing the path in the complex time in two different ways: i) 
First we take a purely imaginary direction, until t = itj such that the trajectory crosses the real manifold Tr = rnK 2 . 
The imaginary part of the angle coordinate of (p^,q^) it then given by —tj/il. ii) Then we start again from (p*,g*) 
and choose the complex phase of each time step dt in such a way that the imaginary part of p remains constant. The 
time t describes then a small loop in the complex time plane that contains the singularity. This gives the order of 
magnitude of the time distance between (jfi,q^) and the singularity, which is in practice extremely small as soon as 
Im(p) is taken reasonably large. For n = 0, there is only one independent (i.e. up to symmetries) singularity, and the 
imaginary part of its time coordinate is just half of t° , the imaginary time required to go from n [— tt, tt] x [—it, tt] 

to Tr n [iT, 3tt] X [— 7T, tt]. 

Such a procedure can be reproduced for various orders n of the integrable Hamiltonian, and we have performed 
it explicitly up to n = 3. Although the method we apply is basically the same, a few important differences may be 
noticed 

i) The number of singularities (i.e. more precisely, of asymptotic regions of the manifold) increases with the order 
of the Hamiltonian. Counting only the number of independent singularities, that is the ones that cannot be 
deduced one from each other by a symmetry, there is only one for n = 0, but (2n + 2) for n — 1, 2, 3. 

ii) If one starts form a point (p^,q^) in an asymptotic region such as Eq. ( [92] ) and propagates along a time path 
that describes a small closed loop of infinitesimal radius around the singularity in time plane, one can show that 
the real part of the resulting momentum is not Re[p*], but Re[p^] + 2-jT/ft, where ft is an integer which depends 
on the order n of the integrable approximation and on the singularity under consideration {ft is equal to one for 
#(°) and to two for and four of the singularities of H^ 2 \ but to three for the two remaining ones). If 
one identifies p and p + 2tt, this means that for n ^ 0, 1, (p, q)(0) are not meromorphic functions. Instead, the 
singularities are of logarithmic type. More precisely, there are ft distinct sheets of the manifold (p, q)(8) around 
each singularity. 

iii) As a consequence, when one tries to reach the complex torus from the neighborhood of a singularity, one should 
specify on what sheet one places oneself. Moreover, this implies that not all singularities are "visible" from the 
real torus: assuming the best way to compute the Fourier integral Eq. (|64| ) is to shift the integration contour 
in the imaginary direction, the only singularities that will be encountered in this way are the ones that can be 
reached by purely imaginary time propagation from the real manifold. For n — 3, only four out of the eight 
singularities are "visible" from the real torus. 

iv) Starting from the neighborhood of a "visible" singularity and following the Hamiltonian flow, one may, depending 
on whether time runs in the positive or negative imaginary direction, and depending also on the chosen sheet of 
the manifold, cross the real manifold Tr in different cell [(21 — l)ir, (21 + 1)tt] x [(21' — X)tt, (21' + 1)tt]. Depending 
on the final cell, the time can be ±ift or ±i(t a — ft) 

On Figs, [l^ and |l3|, we plot, for the resonant torus 10:1 and as a function of the perturbation parameter r, the 
imaginary part of the time coordinate of the "visible" singularities of H^ n > for n — 2 and 3 respectively. What we 
are waiting for is that the k dependence of the V£ :s (for k — rm as well as k 7^ rm) is given by an expression like 
Eq. ww, with t r - s the imaginary part of the time coordinate of the singularity closest to the real torus. On Fi g \L2\ 
and OT are also shown the values 6% obtained by fitting the numerically obtained V£ :s with the expression Eq. Q70]). 
We observe that for n = 2, 3 the variation of the fitted values follows nicely the predicted ones. For higher n, and up 
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to n = 6, the V£ :s are insensitive to the variation of the order of the approximation, and therefore the fitted values 
remains on the curve corresponding to the n = 3 closest singularity. 

The data shown in Fig [l^ and 13 give a pretty convincing picture, which justifies to use confidently the expression 
Eq. ( |70| ) to describe the behavior of the Vf'£ coefficients. Although we believe this to be true from a practical point 
of view, one should, however, resist the temptation to oversimplify this issue. Indeed, it is, to start with, a priori not 
obvious to justify on a rigorous basis the form (mr) 7 V^ :s we have written for the prefactor, and this can only be taken 
as a sensible parameterization. Moreover, even if we did not extend the analysis of the location of all singularities 
for n greater than three because the approach described above becomes somewhat cumbersome, it is still possible to 
locate the closest singularity by a brute force search in the complex 9 plane. Doing this for n = 4, 5, 6 for the the torus 
10:1 at r = 1 shows that this closest singularity slightly drifts as n increase, and that its imaginary part goes from 
0.8 for n = 3 to 0.6 for n — 6, in spite of the VJJm being not affected by this change. This drift, although moderate, 
is still larger than the numerical precision of our fit. In accordance, if we take f(p, q) = cos(p) as done of Fig. [l4|, 
we see that, contrary to SI r - s , the Fourier coefficient of f(8) changes with n even when this latter is greater than 
three, and in particular follows the asymptotic slope 0.6 for n = 6. This indicates that although the basic picture 
we gave to interpret the asymptotic behavior of the is certainly correct, the complete description is presumably 
more complicated and might involve, for instance, the link between the kicked Harper map T and the integrable 
Hamiltonians ff(„) , as well as a more careful analysis of the different ranges in the asymptotic behavior of the V£'£ ■ 
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V. CONCLUSION 



We have presented a scmiclassical picture of how tunneling proceeds in ncar-intcgrable systems. Nonlinear reso- 
nances are identified to provide a key coupling mechanism that leads to a strong enhancement of the tunneling rate 
with respect to the corresponding integrable system. In an analogous way as in the generalized pendulum, these 
resonances induce a periodic modulation of the phase space structures in their vicinity, which couples the invariant 
tori that are symmetrically located with respect to the resonance. Via the succession of one or more such couplings, 
the quantum state may be connected to higher excitations within the well, from which it is significantly easier to 
tunnel outside of the island. 

The practical calculation of tunneling rates between separate wells in near-integrable systems basically includes 
the following four steps. First, the major resonances, with respect to a given value of h, need to be identified within 
the well regions (which at rather large h can be done by direct visual inspection of the Poincare surface of section) . 
Second, the coupling coefficients V£' 8 associated with the effective pendulum Hamiltonian of each of these resonances 
need to be calculated. This requires in general the construction of an integrable Hamiltonian which approximates 
the dynamics with good accuracy. The low-order coupling coefficients are most conveniently evaluated by means of 
the method that was presented in section [II C - i.e., by a Fourier analysis of the action difference SI r:s (9) (Eq. (|64|)) 
between the backward propagation of the resonant torus and the resonant torus itself - whereas high orders coefficients 
can be deduced from the knowledge of their asymptotic behavior. With these coupling coefficients, the transition 
amplitudes across the resonances can be evaluated by means of first, or higher, order quantum perturbation theory 
for rather large, or intermediate, values of h, respectively, and are described by a local semiclassical tunneling process 
deep in the semiclassical regime. As a final step, the transition amplitudes associated with the resonances are then 
combined to form the coupling scheme that the system undergoes to tunnel into the symmetry-related well. For small 
h, many tunneling paths, characterized by the succession of resonances used in the process, are contributing to the 
tunneling rate. As a general rule, we state that the dominant one is such that a maximum number of resonance are 
implied, provided the condition Eq. ( |82| ) is met for all of them. 

To insure that the tunneling mechanism we propose not only is the dominant one, but "makes sense" to start with, 
requires some understanding of the qualitative properties of the coupling coefficients V£ :s . This in addition makes it 
possible to give general criteria concerning which resonances to use in the tunneling path, and how to deal locally 
with the effect of a particular resonance. Such qualitative properties can be obtained from the study of the analytical 
structures, in complex phase space, of the invariant manifold of the integrable Hamiltonian. We have here made a first 
step toward the clarification of these issues, from which some confidence in the global coherence of the description can 
be obtained. We have however not tried to address some more delicate points, such as the computation of prefactors 
in the asymptotic dependence of the V^ s . Moreover, although we have verified for the kicked Harper model the 
independence of the V£ :s on the precise choice of the integrable Hamiltonian, and in particular on the order n of 
the approximation, we can only provide an educated guess on the reason why this property holds. More research 
on complex classical dynamics from a general, mathematical point of view is certainly required to obtain a deeper 
understanding of these issues. 

Despite these remaining questions, we could show that the resonance assisted tunneling mechanism provides a 
quantitative description of tunneling in the near-integrable kicked Harper Hamiltonian. Though exemplified only 
within this particular model system, the theory is described in such a way that it is straightforwardly applicable 
also to other semiclassical one-degree-of-freedom systems with time-periodic perturbations (such as the driven double 
well dynamics [§,[30| for instance). A further generalization that is straightforwardly achievable is the one to decay 
problems - i.e., to the dynamics of quasi-bound states that are coupled to a continuum via a tunneling barrier (as is 
the case e.g. in the a decay of nuclei or, to mention also a more recent example, in laser- or microwave-driven atoms 
P]). Here, the major modification with respect to level splitting problems to which we have restricted our study 
in this paper concerns the spectral observable of tunneling - which would be given by a finite width of a spectral 
resonance instead of a difference of quasidegenerate eigenlevels. In the same way, our theory should apply up to minor 
modifications to the case of autonomous two-degrees-of-freedom systems. 

The quantitative success of resonance-assisted tunneling in the near-integrable kicked Harper model gives us the 
confidence that the underlying mechanism plays a key role also in a more mixed regular-chaotic situation where 
macroscopic parts of the phase space are occupied by chaos. Specifically, we believe that the classically forbidden 
coupling from a regular island into the chaotic sea is now amenable to a natural description by means of the major 
resonances that characterize the dynamics within the island. This still requires to resolve a number of open questions, 
in particular related to the (more or less smooth) transition from regular dynamics inside to chaotic dynamics outside 
the island. We estimate, though, that resonance-assisted tunneling could develop into a powerful tool to quantitatively 
understand semiclassical tunneling processes in mixed systems. 
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APPENDIX A: INTEGRABLE APPROXIMATION FOR RAPIDLY DRIVEN HAMILTONIANS 



The integrable approximation of periodically driven Hamiltonians, such as the the Kicked Harper, is most conve- 
niently constructed by the Lie transformation method. Following its description in standard textbooks |2!j, the Lie 
transformation method introduces a "pseudo" Hamiltonian w e = w e (p, q) in order to perform the canonical trans- 
formation from the original (in general multidimensional) phase space variables (p, q) to the new variables (p, q): 
The evolution from (p,q) to (p, q) is generated by Hamilton's equations of motion associated with w e , where the 
"evolution time" e is given by the perturbation strength of the actual Hamiltonian TL — 7Y(p,q) under consideration. 
The aim is to construct w e in such a way that Ti. acquires, for each value of the perturbation parameter e, a specific 
form (e.g., Tt = T~t{p)) within the new phase space variables (p, q). 

Defining the "inverse evolution operator" T" 1 that undertakes the associated transformation of phase space distri- 
butions according to 

H(p, q) = (IT^Xp, q) := H(p, q) (Al) 
for all (p,q), we can straightforwardly verify that T~ x satisfies 



with the "Lie operator" L f defined by 



J^T" 1 = L e T^ (A2) 



L = dw 1 d__dw ± d_ 
9q dp dp <9q' 



A power series expansion of the original Hamiltonian in the perturbation parameter e, 

oo 

W(p,q) = X> B Wn(p,q), (A4) 
as well as analogous power series expansions of w e , L e , T e , H, 

oo 

w e (p, q) = X] e " w "+i(P' q)i ( A5 ) 

oo 

L e = J2^ L n+i, (A6) 

OO 

T e - 1 = X) enT «~ ) with T o M =Id, (A7) 

n—0 
oo 

H(p,q)^a(p,q) (A8) 

n=0 

can now be used to successively determine H n as well as the coefficients w n of the generating Hamiltonian w e . 

In the case of a one-degree-of-freedom system with rapid periodical driving, we consider a Hamiltonian of the form 

H = H[p, q, 6) = H(p, q,9 + 2tt) with 9 = —t. (A9) 

T 

Here, the small perturbation parameter e essentially corresponds to the driving period t, and the aim of perturbation 
theory consists in e lim inating the explicit time dependence of the Hamiltonian. In order to formally obtain a power 
series of the type (A4), we consider the equivalent, autonomous Hamiltonian system with two degrees of freedom 
given by 

H = H(p,q,J,6)=n (J)+TH 1 (p,q,9) (A10) 

with 



2G 



Ho(J) = J, 
H 1 ( P ,q,6) = ±-H(p,q,< 



(All) 
(A12) 



Here we have additionally introduced the canonically conjugate variables (9, J), corresponding to the scaled time, 
9 = 27rf/r, as well as to the inversely scaled energy, J = —tH/(2it), respectively. The new, time-independent 
Hamiltonian TL can then be written in the form 



H=H(p,q, J) = J+—H(p,q), 

yielding H{p, q) as integrable approximation of the rapidly driven system (| 
Making the power series ansatz 



(A13) 



H{p,q) = J2r n H n (p,q), 



(A14) 



n=0 



we obtain, when applying the Lie transformation method (Al - |A§|) to the Hamiltonian ( A10 ) with p = (J,p), 
q = {9, q) and equating terms of equal power in t, 



n dw n 
2n ~df = " 



n-l 



H n —i — L n - m H m -i — T n _ 1 H 



(A15) 



from the derivative of (Al) with respect to r, and 



m— 1 



n-l 



T(-) = - V L T 



(-) 



with 



m=0 



dwr, d dw„ d 



dp dp dq 



(A16) 



(A17) 



from the relations ( |A2| ) and ( |A3| ), respectively. Requiring w n to be a 2ir periodic function in 9, it follows that we have 
to choose 



u n (p, q,e)=-— IJ2 L n-mH m -i + Tt\H (p, q, 9') d9' 

n J® [m=l ) 



(A18) 



in order to obtain a time-independent expression for H n -\: 



H n -i(p, q) = -(Yl L n . m H m ^ + T^_\h\ (p, q). 



(A19) 



Here we have (according to Q) introduced the notation 

(H)(p,q):=±-[ H(p,q,6)dd, 
{H}(p,q,9) :=H{p,q,6)-{H)(p,q). 

In zeroth order in r, we obtain (with Xjj = Id) 

Wl (p,q,9) = -i- j\H}(p, q ,9')d9' 



(A20) 
(A21) 

(A22) 



yielding the familiar result that the time-periodic Hamiltonian H tends for r — > to its average within one driving 
period: 



27 



H (p,q) = (H)(p,q). 



(A23) 



The expansion coefficients w n , T^_\, H n -i with n > 1 can be successively determined in ascending order in n from 

the expressions ( A16 - A19 ), as the latter depend only on w m , T < ^} 1 , H m -i with 1 < m < n. 
The method shall be exemplified for a Hamiltonian of the form 

H(p,q,9) = H (p,q) + K{6)V{q) 

which typically appears e.g. in the case of electromagnetically driven atoms. Defining 



Kt{9) 



1 

27 



{K}{6')d6\ 



K2{9) :=±- I {KiW)dff, 
27T Jo 

we obtain up to second order in the perturbation r 

dV dH Q 



H{p, q) = H (P, q) + (K)V(q) - t(Kx) 



dq dp 



i (Vi> 2 + \{K 2 ) (K 2 K) 2(K 2 )(K) 
(K 2 



i)\ \ 2 d 2 H 
dp 2 



V dp 



dq \ dpdq dp dp 2 dq 



For the particular case of the Kicked Harper Hamiltonian Eq. ([!]), we have 

oo 

K{0)= S(6/(2n)-n). 



(A24) 

(A25) 
(A26) 

(A27) 



(A28) 



This yields 



T T 

H(p, q) — cosp + cosq — — sinp sing — — (cosp sin 2 q + cosq sin 2 p) 
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sin(2p) sin(2g) + 0(t 4 



(A29) 



as integrable approximation up to third order in r. As evaluated by means of the power series coefficients T„ \ the 
corresponding transformation to the original phase space variables (p, q) reads 



p = p + (t — t) sin q + ~^t{ T ~ i) sin p cos q + ~^i T ~ t) 
x [2(4r — 3t) cosp cosq sing + (t — 2t) sin 2 p sing] , 

- 1 / x - . - 1 / 
q = q - -t(r - t) cosp sing + —t[r - t) 

x [(5r — At) sinp sin 2 q + 2(t — 2t) cosp sinp cos g] 



(A30) 



(A31) 



for < t < t. At times t n — nr — e with integer n and e — > 0+, we have (p, q) = (p, q). This implies that the Kicked 
Harper map Eq. (^) , which essentially corresponds to the stroboscopic section at t = t n of t he dynamics generated by 
H, can be directly approximated by the propagation of the integrable approximation (A2£). 

The terms beyond the third order in r are most conveniently derived with symbolic programs such as MAPLE. The 
coefficients C nm that enter into the corresponding Fourier series expansions C nm cos mp cosng (or sin mp sinng) 

can then be used to propagate the classical equations of motion associated with H as well as to create the quantum 
Hamiltonian matrix to be diagonalized in order to determine the eigenvalues and eigenvectors of H. 
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APPENDIX B: ANALYTIC CONTINUATION OF THE KICKED HARPER TORI INTO THE COMPLEX 

DOMAIN 



In contrast to integrable systems, the analytic continuation of the Kicked Harper tori into the complex domain 
cannot be performed by imaginary time propagation - which, to a certain extent, can be seen from the fact that these 
tori are visualized in phase space by a stroboscopic map, rather than by a smooth Hamiltonian flow. The map itself, 
however, which can be written as 



Po 
Qo 



(po, Qo) 



pa - t sin q 

qo + rsin(p - rsing ) 



(Bl) 



is amenable to complexification, as it represents an analytic vector field in the phase space variables (p, q). In analogy 
to the unbroken KAM tori of the real phase space, the analytic continuations of these tori arc then given by invariant 
manifolds of the complexified map (Bl). This implies that the recursive application of (Bl) on any complex phase 
space point (po,qo) that is located on such a manifold yields a closed, quasiperiodic trajectory which, as in the 
integrable system, is characterized by the same frequency Q (or winding number) as the real KAM torus from which 
on the manifold is departing. 

This fact can be used to construct the analytic continuation of the KAM tori by an optimization program that 
is based on computing complex phase space trajectories with the complexified map (Bl). For this purpose, it is 
convenient to introduce, for each of such trajectories, the "pseudo angle" 



n Q t (mod 2n) 



(B2) 



of the nth iterated point with respect to a given starting point (po, Qo), where f2 is the round-trip frequency for which 
the analytic continuation of the torus is to be constructed. Sorting the trajectory with respect to this pseudo angle 
yields then a smooth, closed orbit - if and only if the trajectory is quasiperiodic with frequency fl. If there is a slight 
mismatch between the frequency of the trajectory and f2 or if the trajectory is not closed at all, then the sorted orbit 
exhibits wild fluctuations, which means that huge differences of (p, q) may occur between adjacent values of 9. A 
phase space point (j>o,Qo) on the "good" manifold may then be found by minimizing a function / = f{po,qo) that 
measures the "noise" of the sorted trajectory that starts from (po, qo) for a given number N of mappings. Specifically, 
we chose / as the average squared difference between adjacent slopes in the sorted orbit - i.e., 



1 

N 



N-l 

E 



Pp{i+1) — Pp{i) Pp(i) — Pp(i-l) 



Mi) 



Mi) 



i=l 

Qp(i+l) ~ Qp(i) Qp{i) - Qp(i-1) 



'p(i-l) 



?p(i+l) 



Mi) 



Mi) 



Mi-1) 



(B3) 



where p{i) represents, for given N , the index permutation that sorts the pseudo angles Qi (i.e., Q v u\ < @pu) for all 
j < i). Clearly, this noise function / will become minimal (and will actually vanish in the limit N — > oo) if all 
phase space points lie on a smooth complex curve - in which case the slope between adjacent points is locally almost 
constant. 

Since / contains many local minima in general, the optimization method has to be implemented in a stepwise way. 
We first start with a r athe r small number N of mappings, corresponding to two or three round-trips of the trajectory. 
Calculating the noise (B3) for a given starting point (p , q ) and for points in its immediate environment yields then, 
by means of the Newton method, a rough (and, at this stage of the algorithm, rather imprecise) estimation of the 
correct phase space point on the KAM manifold. Taking this point as new starting point, the procedure is repeated 
with a higher number N of mappings, which makes the local "valley" around the minimum steeper and thereby allows 
for a more precise estimation of the correct point. The iteration is continued until the point on the manifold is known 
with the desired precision: At each step of the iteration, the starting point (po,Qo) is improved by one step of the 
Newton method, and the value of the noise at the improved point gives an indication by which amount ./V may be 
increased without leaving the "good" minimization valley. At a given maximum number of mappings, convergence of 
the method is checked by comparing the noise f(j>o, Qo) of the optimized starting point (po, qo) with a suitably chosen 
reference value / cr ,t (recall that / should be zero on the KAM manifold for N — > oo). Cases of doubt may be resolved 
by directly looking at the sorted trajectory plotted in pseudo angle space (as shown in Fig. |l5| ). 

The optimization method is additionally simplified for the specific case of the Kicked Harper. Due to the symmetry 
of the phase space with respect to the diagonals (see e.g. Fig. |J), we can infer that the analytic continuations of the 
real KAM tori pass all through the complex hyperplane defined by p = q* . We may therefore restrict our starting 
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points to that hyperplane. By additionally fixing the imaginary part of p (which, in some sense, defines at which 
distance from the real phase space we want to calculate a point on the manifold), only one parameter, namely the 
real part of p and q, needs to be varied for the Newton method. 

A typical example for the analytic continuation of a near-integrable kicked Harper torus was shown in Fig. |^. In the 
corresponding intcgrable system, the two manifolds shown in the figure would join each other across the separatrix to 
form a smooth complex surface. For the near-integrable Kicked Harper, however, the analytic continuations can be 
constructed only up to a certain distance from the real phase space (namely up to Imag(p) ~ 0.8 within the p = q* 
subspace); beyond that value, the optimization program does not converge any longer. A closer look on the outermost 
part of the converged manifold reveals significant fluctuations, as can be seen in the inset of Fig. [jj. Plotting one of 
the last converged trajectories as a function of the pseudo angle (Fig. [l5|) reveals that this trajectory exhibits periodic 
modulations which seem to appear on many different scales. This observation is in fact consistent with the findings of 
Greene and Percival in the standard map, where a fractal structure was found for the outermost part of the complex 



torus (compare in particular Fig. 3 of fll5f with our Fig. 15). It thereby supports our conviction that the reason why 



our method fails to converge is indeed the fact that the natural boundary of the torus is approached. 



30 



APPENDIX C: TUNNELING IN THE PENDULUM 



In this appendix, we give a semiclassical expression for tunneling in the generalized quantum pendulum 

(T — i ) 2 

H = ± + Vi cosrG . (CI) 

2/i 

We shall moreover assume antiperiodic boundary condition for the angle O as a remnant of the existence of Maslov 
indices in the original (p, q) variables. 



1. Invariant torus geometry 

It turns out to be more convenient to adopt an action representation, that is, to work in the Fourier space of the 
pendulum problem. Considering an eigenmode of energy E k , this implies that semiclassically, the manifold Y we are 
interested in is the solution of 

H(I,e) = E k (C2) 



for real action /, but possibly complex angle O. Introducing I+± — Iq± y/2fi(Ek — V\), and I±2 — Jo± y/2[i(Ek + Vi), 
the abscissa of the turning points, on can, as schemed on Fig. [16|, distinguish two parts in this manifold. For / G [/i, ^2] 



for instance, the 2r solutions of Eq. (C2) are real, and we shall denote by the graph of 



9|± 2] (/) = y ± J aiccos [(£ - (/ - / ) 2 /2^)/^i] (i = 0, . . . , r - 1) . (C3) 

This, together with the rj^ 2 _ 1 j defined in the same way in [1—2,1— 1], corresponds to the classically allowed region 

(solid line in Fig. |l6| ). Outside these intervals, has a non zero imaginary part, and the 2r sheets of T come by pairs 
corresponding to complex conjugated functions 

9 [-i,i]( 7 ) = ^±-argch[(£; fc -(/-/ ) 2 /2M)M] (/ = 0, . . . , r - 1) . (C4) 

for Ie [I-i,h], 



Qi [2 oo] GO = ~ — ± -argch [(E -(I- Io) 2 /2^)/^i] (I = 0, . . . , r - 1) . (C5) 

for / > I2, and ©ji^ _ 2 ]C0 defined in the same way for / < 1-2- The interrupted line in Fig. |l^ represent a 
complex path ©(/) (I real) in this classically forbidden region. In the following, we will denote by r|^ b j ([a, b] = 
[-co, —2], [—2, —1], [—1, 1], [1, 2], [2, 00]) the various parts of T described by the ©J^ above. We shall also make use 
of the notation l[ a ,b](0) for the inverse of ©J~ b ](-0 on the interval € [(22 — l)ir/r, 2ln/r[ and of ©{^yCO on the 
interval 6 e [21-ir/r, (21 + l)%/r\. 



2. Local solutions and connection formula 

Following standard WKB analysis, a local approximation to the Schrodinger equation H 1 ^ — E^ can be obtained 
between two turning points, ie. under each of the r|^ fc j as 

^M]( J ) = /-,+ ,, exp ?r — ' (C6) 



with the action integral 



U i el ^ {I) rJI laM 



% = e^(6)de, (C7) 
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( 5[ofc](^) = |7fcVirsin ( r ®{^6](-0J I' an< ^ the period of the classical motion at energy E^. 

This analysis can be legitimated a priori when the interval spanned by the manifold on the I coordinate is large 
as compared to the typical domain where the WKB construction diverges at the caustics. Here it means that 
|7a — ii| ^> h 2 / 3 , which can be written 

# » 3 (C8) 

if V\ is small as compared to energy Ek ■ 

The first step in the construction of a global wavefunction is then to patch together these local solutions at caustics. 
A convenient way to do this here is to use Langer |38| connection formula which, for a given choice of the origin of 
the action integrals, imposes some relations between the various coefficients A. Let us denote for instance by A Q ut 
and Am the amplitude of the outgoing and incoming wave on the classically authorized side of a turning point, and 
by A > and A < the ones of the exponentially decreasing and increasing functions on the classically forbidden side. 
Then, except for a global multiplicative constant A, the relation between these various amplitudes can be expressed 
in term of a single parameter ry. Assuming the origin of the phase integral is taken at the turning point this relation 
can be expressed as 



Aut = -4expH(7r/4 - rj)} (C9) 

A„ = -4exp[+i(7r/4 - ??)] (CIO) 

A> = AcQs[r}} (Cll) 

A < =A2sm[r]}. (C12) 



Neglecting tunneling amount to assume that only the exponentially decreasing function should be used, in which case 
the dephasing between the incoming and outcoming wave on the classically allowed side of the caustic is tt/2. The 
above expressions generalize this to the case where there is some non zero amplitude on the exponentially increasing 
wave. 

For the pendulum the parameter r\ is obviously zero for caustics at I = I±2 since there, any amplitude on an 
exponentially increasing function would lead to a non normalisable solution. Moreover, by symmetry, the r caustics 
at 7 = 7i, (respectively J_i), must have the same parameter r\ = r\\ (respectively T)-i)- 



3. Paths within the classically allowed regions 

To obtain a global approximation to the Schrodinger equation, it is now necessary to choose the manifold T, that 
is the energy Ek, and the caustics parameter 771 and 7j—i, in such a way that the various connections are consistent 



one with each other, or in other words that, once a choice of the origins of the action integrals in Eq. (C7) is made, 
a given set of parameters A 1 ^ fulfill all the constraints due to the connection formula. This is equivalent to the fact 
that such condition apply to any independent loop of T. 

The first closed path of the manifold that we are going to consider is the one belonging to the classically allowed 
region [Ti,^]- Moving in the upward direction on this loop, one gains a phase —n/2 + 7/1 at each of the r turning 
points with abscissa 7i, and +7r/2 at each of the one at I = 7a . Consistency should therefore imposes that 

Stat I h + 2r77i = tt [mod 2tt] . (C13) 

where S^f is defined as in Eq. (|C7|), but between (J, 6) and (I, 6 + 2tt). 

One should be a little bit more careful here because the geometry of the phase space of the pendulum is a torus, 
— [121 

which implies that S\. ot is not independent of the origin (7, 9) of the integration. Indeed, the fact that the final point 
is really (7, 9 + 2tt), and not (7, 9), implies that dS [ ^/dI = -2tt, and that therefore s[^ ] (I) = s[^ ] (7 ) - 2tt(7 - 7 ). 



It is therefore impossible to impose the condition Eq. (|C13| ) for all 7. However, because the phase space is a cylinder, 
the only admissible values of the action are, assuming antiperiodicity in 9, the discrete set I n = h(n + 5), with n a 



positive or negative integer. For all of those, Eq. (C13) remains unchanged if we replace Stot by <Stot = Stot + 27r7 
(since 2irl n /h is equal to zero modulo 27r for all 7„). One therefore recovers in this way the "natural" quantization 
condition 

S l t H ] /h + 2r Vl = jl 7 [12] (9)d9 + 2r m =^(k + ^\ . (C14) 
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Using the same reasoning, we obtain also that 



S [ t - t 1 ~ 2] /h-2rr ] _ 1 = ^I [ _ 1 _ 2] (e)dQ - 2rr,-i = 2tt (k> + \ 



(C15) 



4. Classically forbidden paths 

The second kind of closed loop that we shall consider is, as usual in tunneling problems, the one starting from some 
turning point (Ji, 0°), with 0° = for some given I = 0, . . .r — 1, following one of the complex branch of T up to 
(I-i, 0°), switch to the complex conjugated branch and move back to (2i, 0"). Consistency then impose that 

tan?7i tan?7_i = ^ cxp (— 2^\ , (C16) 



4 * V h 

with 

rh i r Ii 



* = L\ |Im (e[ - Ml(/)) W = \ I [ argch (t - rf/ ■ (C17) 

Finally the last kind of loop that we need to consider is the one which starts from some turning point Qf), follows 
one of the complex branch of T up to (i-i, 0°), moves upward on the classically allowed branch up to 0?+i), 
then backward on a complex branch to Q^+i), and finally down to (ii, 0°). Assuming Eq. ( C16| ) is verified, which 



make unimportant which of the two possible complex branches is used at each traversal of the classically forbidden 
region, this still imposes a constraint on the phase of the wavefunction, namely, after an integration by parts of the 
action integral 



1 f °+i 



• / 1+1 dO [I [li2] (0) - J M ,_ 2] (0)] -2(r 7l +r?_ 1 ) =0 [mod 2tt] . (C18) 
Je? 

Any other loop on T can be decomposed onto the ones considered above. 



5. Quantization conditions 



The left hand side of Eq. (C16) being exponentially small in h, this implies that either r]i or r/_i are exponentially 
small, but, assuming 



/, i = // ( K + - ) ■ 



(C19) 



where kq is not an integer, not both. We assume for instance that 771 -C T)-i ~ 1. 

This means that the real torus T[ 12 ] almost quantizes in the absence of the corrections due to tunneling, ie. that 



sl 12] /h=^i [12] (e)de^2n {k+ l - 



(C20) 



The normalization of the wavefunction furthermore implies that the modulus of all the -4f 12 i is one. 

Our goal her e is n ot to e valua te the shift in energy associated with the modification of the quantization condition 
between Eqs. (C14) and ( |C20 ). Indeed this is presumably much smaller than polynomial corrections in H that we 



have neglected. Rather, what we are interested in are the amplitudes A|^ 2 _i] °f the part of the wave function which 



has tunneled on the other side of the resonance. We thus assume that the energy fulfills Eqs. 
that the manifold T is fixed. Neglecting again 771 , the parameter ry_i is then determined by Eq. 





h J e o +1 de [ J[i < 21 (e) _ /[ - x '- 2] (e) ] [mod 2n] 



^20D , and therefore 
as 

(C21) 
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Note that this choic e of ?7_i necessarily fulfills Eq. ( C15| ), but that the reciprocal is not true. Indeed, as we shall 
see below, Eq. ( |C18|) imposes the "selection rule", which stands that only W.K. B. modes whose quantum numbers 
differ in a multiple of r are coupled. Using the quantization condition Eq. (C14) and that I[i,2] (O) and Jr_i ) _2](@) 
are symmetric with respect to Iq — h(no + we obtain 



V-i 



2ir(k — kq) 



(C22) 



Applying then Eq. ( C16| ) fixes the value of rji as 

1 



4tan?y_ 



exp -2- = 



4 tan 



2ir(k — k,q) 



exp I —2 



(C23) 



and, taking the origin of the phase integral at 0) for all the S l ^ 2 _ 1 i (and at 0) for the Sfj^i) the transmitted 
amplitude as 



A 



l t 2 -i\ - a t exp (-i(2l ± l)^ ± wr/4) 
1 

At = , . r^r- exp 



. / 2ir(k - k ) 
2 sin - 



(C24) 
(C25) 



6. Selection rules 



Going now back to the representation, and evaluating the corresponding Fourier transformation at the stationary 
phase approximation, we obtain 



*(0) 



1 



exp -5 [1)2] (6) + At exp -S[_ 2 ,_i] (6) - i(2z/(6) + 1)tj_i 



(C26) 



with v(Q) = Int[0/(27r/r)]. It is straightforward to generalize this expression to a more general potential, like in Eq. 
( |30| ) , as long as the general structure of branches, the periodicity and the symmetry are preserved. 
Now we introduce the quasimodes 



u k >(e) 



T k <\G k 



:exp(^(6)/fr) 



(C27) 



which are built on real invariant tori k' such that 



S t k f = JJl kl {Q)d@ = 2'K (* + \ 



(C28) 



'J is essentially equal to u k , plus some exponentially small tunneling contribution ^t localized on the other side 
of the resonance. If we write the scalar product of ^t with some u k >, assuming — (k' — kq) — {k — kq) (and thus 
T k , ~T k = T, and \Q k \ ~ |6 fe /| = 9) we obtain 



( Uk , |* T > = ^ J ^ exp i (5 [ _ 2 ,_ 1] (Q)/h - (2/(6) + l)?/-i - ft, (Q)/h) 

r r 27T/r de 



A T x 



TJ 



e 



expi (^^(GJ/ft-Tj-i - Sfc/(0)/ft) 



1 r " 

- ^ exp (S[_ 2 ,_i] (27r/r)/fi - 2r?_! - ft. (27r/r)/fi) 



(C29) 



(C30) 



34 



The integral on O can be evaluated using that S[—2,—i\ (@) — Sk' (©) — t(@ )(Ek — Ek>)- For the remaining sum on /, one 
should note that Sk'(2Tr/r)/h — 2itk' /r, and that in the same way, Eq. ( |C18| ) imposes that iS , r_ 2i _i](27r/r) — 2t?_i = 
S'[i2](27r/r) [mod 2n] = 2kn/r [mod 2tt\. We therefore obtain that, up to a global phase 

<«* i*t) = ^ sin( f F f; A) e . ( csi ) 

TroE/rA 

' m 

with i5_E = — and A = /iu; is the spacing between the quasimodes energy at E ~ Ek — E' k . Eq. ( C3l| ) contains 
the selection rule, stated above, that a resonance with r islands couples through tunneling quasimodes with quantum 
numbers differing in a integer multiple of r. 



7. Comparison with the result of perturbation theory 

The semiclassical treatment described can be justified a priori in the domain ft 2 / 3 <C ^fT^^y/Ek + V\ — \/Ek — Vi), 
with Ek — h 2 (k + i) 2 /(2/i) which is disjoint from the regime of validity of the quantum theory of perturbation. 
However we can show that for 

V 1 « (/ fe -2o) 2 
\h~h\ > fir (C32) 
?7_i < l 1 . 

the two formulas coincide. Indeed in this case the transmitted WKB function has an overlap with the plane wave 
|fc — rm) almost equal to 1, and then we show in the following that A T ~ (fc — rm^" 1 - 1 ), with the estimation of 
the k th eigenmode of the system given by the theory of perturbation at the order to. Whether the coincidence holds 
for every r\ is not obvious for now. 

The m th order perturbation theory gives 



/ T , x to to 

^ ' 1=1 



{k — kq) 2 — (k — kq — lr) 2 ' 



(C33) 



">V- -[^Y " K ,1, ..7 = (A)" (C34) 



which, in the limit of large to, is equivalent to 

'2(k-K ) 



h 2 r 2 ) m , r p(fc-K ) ^ \h 2 r 2 ) m!r ( TO+ !Li 
r 



Now we use the Stirling formula and the equivalent xT(x) ~ 1 for x — > 0, and we get 



(*-H* W )^fS) ■ (C35) 



2?7_i \fi, 2 r 2 TO 2 

On the other hand, the classical action a along the imaginary path relating the k th K.A.M. torus to its symmetric 
is equivalent to 



in the perturbative regime. This can be rewritten, by using (C22) and ?7_i <C 1 



exp 

and the expression for At is 



h ( 7r ) \ ^\h 2 r 2 m 2 ) \h 2 r 2 m 2 )) ^ ^ 

Then 

' (C38) 



2r/_! \h 2 r 2 m 2 



which is identical to (C35). 
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FIG. 1. Eigenphase splitting for the kicked Harper system in the nearly integrable regime (r = 1). The splittings correspond 
to the Floquet mode whose WKB approximation quantizes on the classical torus with action 7„/ 4 = 7r/4, for different values of 
N = 2n/h. Squares: exact numerical calculation. Filled circles: semiclassical approximation based on the resonance assisted 
tunneling mechanism. The doted line corresponds to the splitting for the integrable approximation. The solid and dashed lines 
are just guides to the eye. 
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FIG. 2. Poincare section for the kicked Harper map Eq. (|2j) with r = 1. The eigenphase splitting shown on Fig. |l] correspond 
to the tunneling processes between the quasimodes constructed on the tori with action I = pi/4 marked as thick lines in the 
left and right wells. 
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q 

FIG. 3. Poincare section for the integrable approximation with n — 6 of the the kicked Harper map, see Eq. ([Tij) at 

t = 1. Note the good agreement between the invariant tori of this integrable system and the K.A.M. tori of the kicked Harper 
at t = 1 (Fig. |). 
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FIG. 4. Part of the analytic continuation of an invariant torus and its symmetrically shifted counterpart for the integrable 
Harper hamiltonian. The imaginary part of momentum of the manifold is plotted versus the real parts of position and 
momentum, on top of the real phase space. 
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FIG. 5. Part of the analytic continuation of a KAM torus and its symmetrically shifted counterpart in the near-integrable 
Kicked Harper at r = 1. The imaginary part of momentum of the manifold is plotted versus the real parts of position and 
momentum, on top of the real phase space. In contrast to integrable dynamics, the two manifolds on each side of the separatrix 
do not join each other, but are interrupted far inside the regular phase space domain. Note the fluctuations that appear at the 
outermost end of the converged manifold (as shown in the inset for the manifold on the right hand side). 
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FIG. 6. Poincare section in the vicinity of the 10:1 resonance of the kicked Harper, in the action angle variables of the 
integrable approximation H^ n ' (n = 6). Note the similarity of the phase space structure with the one of a generalized 
mathematical pendulum, with 10 islands instead of one. 
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FIG. 7. Visual representation of resonance-assisted tunneling: Plotted is, for the Kicked Harper system with r = 1, the 
succession of complex paths that are used to accomplish the tunneling process for the 10th excited state at 2n/h = 126. 
Instead of a direct transition (dashed curve) to the symmetry related torus, the quasi-mode is coupled via the 8:1 (hardly 
visible) the 10:1, and the 14:1 resonance to the 50th excited state, which lies sufficiently close to the separatrix to tunnel with 
a large rate to its counterpart. 
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FIG. 8. For the tenth iterate of a set of randomly distributed points in phase space, distance, as a function of the order of the 
approximation n, between the sets of resulting points depending on whether the kicked Harper or its integrable approximation 
is used. The perturbation parameter is r = 1.0. 
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FIG. 9. fc-dependence of the resonance coefficients Vj. 10:1 for the resonance 10:1 of the kicked Harper with r = 1. The 
various symbols correspond to different orders n of the integrable approximation used for the calculation, demonstrating 

the independence of the V£ :s on n for n ~ no = 6. The large circle symbols at k = 10, 20, 30 show the result obtained for these 
Fourier coefficients with the pseudo-separatrix method used in Ejl. 
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FIG. 10. Top : same as Fig. [j], but for the Floquet mode whose WKB approximation quantizes on the classical torus with 
action = n/6. Bottom : the classical action, for each value of N, of the modes involved in the dominant tunneling path. 
The horizontal doted line corresponds to the initial action The horizontal solid lines correspond to the actions of the 

resonances 8:1, 10:1 and 14:1. The curved solid lines are such that / — I^/g = const, x h (that is, the distance from is 
Constance in mode number unit). Notice that as N = 2n/h increases, more resonances are involved. 
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FIG. 11. Eigenphase splitting for the kicked Harper system in the nearly integrable regime (r = 1), for a fixed value of h 
(= 27r/100), and a variable initial torus. Squares: exact numerical calculations. Dots: integrable semiclassical approximation, 
without including the effect of the resonances. As the initial torus get close enough from the separatrix J m ax = 7T, the 
integrable approximation for the tunneling reproduces faithfully the exact quantum result, showing that the chaotic layer near 
the separatrix plays no role here. 
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FIG. 12. For the order n — 2 integrable approximation H^ n >: comparison, as a function of the perturbation parameter r, of 
the fitted value of the slope of the log of the coefficient V^ 0:1 (squares) with the imaginary part of the angle coordinate of 
the visible singularities on the torus 7io:i (circles) The lines are just guides to the eye. 
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FIG. 13. Same as Fig. but for the order n = 3 integrable approximation H 
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o Q f(8) = cos p(0) - order 3 
A A f(0) = cos p(9) - order 6 



-10 



-20 




-30 







10 



20 

k 



30 



40 



FIG. 14. Fourier coefficients ft of the functions f(6) = cos [p(8, iio:i)]- The different curves correspond to different orders 
n of the integrable approximation used. We observe that although the fk's remain insensitive to the choice of the integrable 
approximation for not too large k, their asymptotic behavior differ depending on whether the integrable approximation with 
n — 3 or n = 6 is chosen. 
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FIG. 15. Sorted complex orbit belonging to the outermost part of the continuation manifold shown in fig. JE). Plotted is the 
real part of momentum as a function of the pseudo angle 9. One clearly recognizes significant periodic modulations of the orbit 
which appear on different scales, as can be seen from the magnifications shown in the insets. 
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